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Abstract 

Wc observe that the exactly solved eight-vertex solid-on-solid model contains an hitherto 
unnoticed arbitrary field parameter, similar to the horizontal field in the six-vertex model. 
The parameter is required to describe a continuous spectrum of the unrestricted solid-on- 
solid model, which has an infinite-dimensional space of states even for a finite lattice. The 
introduction of the continuous field parameter allows us to completely review the theory of 
functional relations in the eight-vertex/SOS-modcl from a uniform analytic point of view. 
We also present a number of analytic and numerical techniques for the analysis of the Bethe 
Ansatz equations. It turns out that different solutions of these equations can be obtained from 
each other by analytic continuation. In particular, for small lattices we explicitly demonstrate 
that the largest and smallest eigenvalues of the transfer matrix of the eight- vertex model are 
just different branches of the same multivalued function of the field parameter. 
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1 Introduction 



The powerful analytic and algebraic techniques discovered by Baxter in his pioneering papers [1- 
4] on the exact solution of the eight-vertex lattice model laid the foundation for many important 
applications in the theory of integrable systems of statistical mechanics and quantum field theory. 

This paper concerns one of these techniques — the method of functional relations. Over 
the last three decades, since Baxter's original works [1-4], this method has been substantially 
developed and applied to a large number of various solvable models. However, the status of this 
method in the eight-vertex model itself with an account of all subsequent developments has not 
been recently reviewed. This paper is intended to (partially) fill this gap. Here we will adopt an 
analytic approach exploiting the existence of an (hitherto unnoticed) continuous field parameter 
in the solvable eight-vertex solid-on-solid model of ref. [3]. 

For the purpose of the following discussion it will be useful to first summarize the key 
results of [1—4]. Here we will use essentially the same notations as those in [1]. Consider the 
homogeneous eight-vertex (8V) model on a square lattice of N columns, with periodic boundary 
conditions. The model contains three arbitrary parameters u, rj and q = e J7rr , Imr > 0, which 
enter the parametrization of the Boltzmann weights (the parameter q enters as the nome for the 
elliptic theta- functions). The parameters r\ and q are considered as constants and the spectral 
parameter u as a complex variable. We assume that the parameter rj is real and positive, 
< rj < it/2, which corresponds to the disordered regime [5] of the model. 

The row-to-row transfer matrix of the model, T(u), possesses remarkable analytic properties. 
Any of its eigenvalues, T(u), is both (i) an entire function of the variable u, and (ii) satisfies 
Baxter's famous functional equation, 

T(«) Q{u) = f(u - rj) Q(u + 2rj) + f(« + rj) Q(u - 2rj), (1.1) 

where Q 

f(u) = (M^\0i)f , (1.2) 
and Q(u) is an entire quasi-periodic function of u, such that 

Q(u + tt)=±(-1) JV/2 Q(u), Q(n + 2vrr) = q~ 2N e~ 2iuN Q(u). (1.3) 

These analytic properties completely determine all eigenvalues of the transfer matrix T(u). In- 
deed, Eq. ljl.ip implies that the zeroes u\,U2, • • • , u n , of Q(tt) satisfy the Bethe Ansatz equations, 

Ti \ = ~rw o~T' QK =0, k = l,...,n. 1.4) 

f(u k - rj) Q(u fc - 2r/) 

These equations, together with the periodicity relations (|1.3j) . define the entire function Q(u) 
(there will be many solutions corresponding to different eigenvectors). Once Q(u) is known the 
eigenvalue T(it) is evaluated from (jl.ll) . 

The entire functions Q(k) appearing in (jl.ip are, in fact, eigenvalues of another matrix, 
Q(tt), called the Q-matrix. Originally it was constructed [1] in terms of some special transfer 
matrices. A different, but related, construction of the Q-matrix was given in [2] and later on 
used in the book [5]. An alternative approach to the 8V-model was developed in [3,4] where 
Baxter invented the "eight-vertex" solid-on-solid (SOS) model and solved it exactly by means 



1 Here we use the standard theta-functions [6], #»(u | q), i = 1, . . . , 4, q 
7T and 7tt. Our spectral parameter u is shifted with respect to that in [1] by a half of the imaginary period, see 
Sect l2.5.Tl for further details. 
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of the co-ordinate Bethe Ansatz. This approach provided another derivation of the same result 
(|l.ip - (|1.4p . since the 8V-model is embedded within the SOS-model. 

Baxter's Q-matrix (or the Q-operator) possesses various exceptional properties and plays an 
important role in many aspects of the theory of integrable systems. A complete theory of the Q- 
operator in the 8V-model is not yet developed. However for the simpler models related with the 
quantum affine algebra U q (sl2) (where the fundamental L-operators [7] are intertwined by the R- 
matrix of the six- vertex model) the properties of the Q-operator are very well understood [8]. In 
this case the Q-operators (actually, there are two different Q-operators, Q+ and Q~) are defined 
as traces of certain monodromy matrices associated with infinite-dimensional representations of 
the so-called ^-oscillator algebra. The main algebraic properties of the Q-operators can be 
concisely expressed by a single factorization relation 

T+(«) = Q+(n + (2j + 1)^)Q_(« - (2j + l)r?) (1.5) 

where T^(u) is the transfer matrix associated with the infinite-dimensional highest weight rep- 
resentation of U q (sl2) with an arbitrary (complex) weight 2j. Remarkably, this relation alone 
leads to a simple derivation of all functional relations involving various "fusion" transfer ma- 
trices and Q-operators [8,9]. For this reason Eq. (|1.5|) can be regarded as a fundamental fusion 
relation: once it is derived, no further algebraic work is required. 

An important part of the theory of the Q-operators belongs to their analytic properties with 
respect to a certain parameter, which we call here the "field parameter". In the context of 
conformal field theory (considered in [8,9]) this is the "vacuum parameter", which determines 
the Virasoro highest weight A; in the six-vertex model it corresponds to the horizontal field. In 
fact, the very existence of two different solutions [8, 10] of the TQ-equation (jl.lj) can be simply 
illustrated by the fact that the spectrum of the transfer matrix does not depend on the sign of 
the field, whereas the spectrum of the Q-operator does. 

It is well known that it is impossible to introduce an arbitrary field parameter into the "zero- 
field" or "symmetric" eight vertex model of [1] without destroying its integrability. However, such 
parameter is intrinsically present in the solvable SOS-model. It does not enter the Boltzmann 
weights, but arises from a proper definition of the space of states of the model. To realize this 
recall that the SOS-model [4] is an interaction-round-a-face model where the face variables ti 
(called the heights) take arbitrary integer values — oo < l{ < +oo. Its transfer matrix acts in 
an infinite-dimensional space of states even for a finite lattice. It has a continuous spectrum, 
parameterized by the eigenvalue of the operator which simultaneously increments all height 
variables, ti — > ti + 1, on the lattice. Indeed, taking into account the results of [11, 12], it is 
not difficult to conclude that the calculations of [4] require only a very simple modification to 
deduce that the eigenvalues of the SOS transfer matrix enjoy the same TQ-equation (jl.ip . but 
require different periodicity properties 

Q ± (u + tt) =e ± ^Q±(u), Q±(u + 2vrr) = q~ 2N e** e~ 2luN Q±(u), (1.6) 

where the exponent jp is arbitrary. It is determined by the eigenvalue u = e 2triv ^ n of the height 
translation operatocl (the second exponent tp is dependent on tp). It is natural to assume, that 

2 In [4] Baxter restricted the parameter r\ to the "rational" values Lr\ — min + m,2TTT, L,mi,rri2 G Z and 
considered a finite-dimensional subspace of the whole space of states, regarding the values of heights modulo L. 
In this case the phase factors ui = e 2lw '' r take quantized values uj l = 1 (see [11, 13] for further discussion of 
this point). Apart from providing the conceptual advantage of a finite-dimensional space of states, the above 
restriction on r\ and uj was not used anywhere else in [4] and, therefore, can be removed. The transfer matrix of 
the 8V-vertex model (reformulated as the SOS-model) acts only in the finite-dimensional subspace of the SOS 
space of states, corresponding to a discrete set of exponents tp — kn and ip = (the value of N is assumed to be 
even) . 
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the functions Q±(u) solving these equations, are eigenvalues of the Q-operators for the SOS- 
model. Of course, it would be very desirable to obtain their explicit definition (and generalize 
the algebraic result (|1.5j) to the SOS-model), however, many properties of these operators can 
already be deduced from the information about their eigenvalues. 

In this paper we will develop the analytic theory of the functional relations for the SOS-model 
starting from the eigenvalue equations (jl.ip and (|1.6j) . Bearing in mind that the TQ-equation 
arises from very non-trivial algebraic "fusion" relations [1], it is not surprising that it 
implies all other functional relations. The required calculations are essentially the same as those 
in [8,9], apart from trivial modifications arising in the context of lattice models. 

The eigenvalues Q±(u) are two linear independent "Bloch wave" solutions [8,10] of the finite 
difference equation (jl.ip for the same T(-u). Their quantum Wronskian \N((p), defined as, 

2i\N(<p)f(u) =Q + (u + ri)Q_(u-ri)-Q + (u- V )Q_(u + i ] ) , (1.7) 

is a complicated function of <p, r\ and q, depending on the eigenvalue T(ii). The Bloch solutions 
Q±(u) are well defined provided the exponent p does not take some "singular values" (see 
Eq. (|2.17p below), where W(c/?) vanishes. Otherwise Eq. (jl.ip has only one quasi-periodic solution, 
while the second linear independent solution does not possess any simple periodicity properties. 

All singular cases (in fact, they split into different classes) can be effectively studied with 
a limiting procedure starting from a non-singular value of p. In the simplest case, when n is 
generic and ip approaches the points <p = kir, k G Z, the solutions Q+(«) and Q_(u) smoothly 
approach the same value (which for even N coincides with the eigenvalue Q(u) of the 8V-model). 

A more complicated situation occurs when the field <p tends to a singular value, say p = 0, 
while the parameter ij simultaneously approaches some rational fraction of n, where the transfer 
matrix of the 8V-model has degenerate eigenvalues. The limiting value of T(it) is always uniquely 
defined. However, if T(u) is a degenerate eigenvalue, the limiting values of Q±(u) are not 
uniquely defined. They have "complete exact strings" of zeroes whose position can be made 
arbitrary by changing the direction of the two-parameter (77, (/?)-limit. Obviously, this reflects 
a non-uniqueness of eigenvectors for degenerate states [13]. An immediate consequence of this 
phenomenon is that, for rational r], there is no unique algebraic definition of the Q-operator in 
the symmetric 8V- model. This explains an important observation of [14], that Baxter's two Q- 
operators, constructed in [1] and [2], are actually different operators, with different eigenvalues 
for degenerate eigenstates. 

Further, the eigenvalues Q±(u), considered as functions of <p, have rather complicated ana- 
lytic properties. Besides having the (relatively simple) singular points discussed above, they are 
multivalued functions with algebraic branching points in the complex </>plane. The correspond- 
ing multi-sheeted Riemann surface appears to be extremely complicated; we were only able to 
numerically explore it for some particular eigenvalues. 

It is easy to see that the roots of the Bethe Ansatz equations (jl.4p , considered as functions of 
p satisfy a system of the first order ordinary differential equations, duk/dip = Uk{ui,U2, ■ ■ ■ , u n ), 
where U], are meromorphic functions of their arguments. Using these equations one can ana- 
lytically continue any particular solution of (|l,4p along a continuous path between two points, 
corresponding to the same value of p on different sheets of the Riemann surface. In general, the 
resulting set of roots u\,U2, . . . ,u n differs from the initial one, but, of course, satisfies exactly 
the same equation (ll.4p . In other words, different solutions of the Bethe Ansatz equations and, 
therefore, different eigenvalues of the transfer matrix can be obtained from each other by the 
analytic continuation in the parameter p>. 

Guided by the above observation one might be tempted to suggest that the Bethe Ansatz 
equations have only one solution, considered as a function of p. Undoubtedly, that could be an 
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elegant resolution to the problem of completeness, which traditionally attracts a lot of attention 
in the literature (see, e.g., the recent papers [13,15] and references therein). At the moment we 
cannot prove or disprove the above assertion. The analytic structure of the eigenvalues of the 
eight-vertex SOS model is quite complicated and certainly deserves further detailed studies. 

This material was planned as an introductory part for an extended version of our previous 
work [16, 17] devoted to the connection of the 8V-model with the Painleve transcendents. How- 
ever, in the course of writing, we realized that a review of the theory of the functional relations 
in the 8V-model could be of interest to a much wider audience than originally intended and 
deserves a separate publication. A detailed account of the results presented in [16, 17] will be 
given in the second paper of this series [18], which is totally devoted to the special rj = tt/3 
case of the eight vertex model with an odd number of sites. There we will consider remarkable 
connections of this special model with various differential equations, including the celebrated 
Lame, Mathieu, Painleve III and Painleve VI equations. 

The organization of the paper is as follows. In Sect. [2] we present the analytic theory of 
the functional relations in the 8V/SOS-model. In Sect. [3] we discuss some applications of the 
quantum Wronskian relation. In particular, we show how it can be used for the analysis of 
the degenerate states. In Sect. 0] we completely classify eigenvalues of the transfer matrix of 
the symmetric 8V-model for small lattices of the size N < 4. We then study the analytic 
properties of the eigenvalues with respect to the field variable ip with a combination of analytic 
and numerical techniques. In Conclusion we briefly summarize obtained results. Basic properties 
of the 8V-model are reviewed in the Appendix. 

2 Functional relations in the eight-vertex SOS model. 
2.1 Overview 

In this section we will outline the analytic theory of the functional relations in the SOS-model 
(which also covers the symmetric 8V- model) . Actually most of the functional relations discussed 
below are quite universal and apply to a wider class of related model. They include the six-vertex 
model in a field [19], the restricted solid-on-solid (RSOS) model [20] and some integrable models 
of quantum field theory: the c < 1 conformal field theory [8] and the massive sine-Gordon model 
in a finite volume [21]. 

Let T(n) and Q(u) denote the eigenvalues of the transfer matrix and the Q-operator respec- 
tively and 77 is an arbitrary real parameter in the range 



In all subsequent derivations we will use only one general assumption about the properties of 
the eigenvalues: 

We assume that T(u) and f(u) are entire periodic function of the variable u, 



< 7] < vr/2 . 



(2.1) 



J( u + vr) = J(u) 



f( n + vr) = f(u) 



(2.2) 



and that the function Q(u) solving the TQ-equation, 



J{u) Q(tt) = f (it - rj) Q{u + 2r?) + f (u + rj) Q(u - 2rj) 



(2.3) 



is a also an entire (but not necessarily periodic) function of u. 
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For every particular model the above requirements are supplemented by additional, model- 
specific analyticity properties of Q(u) (such as, for example, the imaginary period relation fjl .3j) 
for the 8V-model). These properties will only be used in Sect. [3]and[4j they are discussed at 
the end of this section. 

As explained in the Introduction, once the additional analyticity properties are fixed, the 
functional equation (|2.3[) completely determines all eigenvalues T(u) and Q(u). For certain 
applications, however, it is more convenient to use other functional equations in addition to (or 
instead of) (12. 3D . We will show that all such additional functional relations in the SOS-model 
(and in the related models mentioned above) follow elementary from two ingredients: 

(i) the TQ-equation itself (Eqs. (|2.2I) and (|2.3I) above), and 

(ii) the fact that for the same eigenvalue T(u) this equation has two different [8, 10] linearly 
independent solutions for Q(u) which are entire functions of u. 



The only property of the function f (it) essentially used in this Section is its periodicity (|2.2j) . 
For technical reasons we will also assume that f(u — rf) and f(n + r/) do not have common zeroes. 
This is a very mild assumption, excluding rather exotic row-inhomogeneous models, which are 
beyond the scope of this paper. 

2.2 General functional relations. 

Since T(u) is an entire function, Eq. (|2.3p implies that the zeroes u\, ■ ■ ■ , u n of any eigenvalue 
Q(u) satisfy the same set of the Bethe Ansatz equations 

Hvk + rj) Q(u k + 2 V ) 

?K^) = -QK^)' QK) = < fc = 1 '-.". ^ 

where the number of zeroes, n, is determined by the model-specific analyticity properties. 

For any given eigenvalue T(u) introduce an infinite set of functions T k (u), k = 3,4, ... oo, 
defined by the recurrence relation 

Tfc(ii + rj)J k (u - 77) = f(u + kr])f(u - kr]) + T fc _i(u) J k+1 (u), k>2, (2.5) 

where 

T («) = 0, Ti(u) = f(u), J 2 (u) = T(u) . (2.6) 
This relation can be equivalently rewritten as 

T(u) J k {u + krj) = f(u - 77) T fc _i(n + (k + l)r?) + f(« + rj) J k+ i{u + (k — l)r/) , (2.7a) 

or as 

T(u) J k (u - krj) = f(u + t?) T fc _i(u - (k + l)??) + f(u - 77) T fc+ i(n — (k — 1)? ? ) . (2.7b) 
Using the definition (j 2 . 5 j) one can easily express T k (u) in terms of T(u) as a determinant 

T fc («) = (f (fe) («)) _1 det M ab (u + kri) , k>2, (2.8) 

l<a,b<k—l 

where the (k — 1) by (k — 1) matrix M(u) a b, 1 < a, b < k — 1, is given by 

M ab (u) = 5 a)b J(u - 2arj) - 5 aM1 f(u - (2a + 1)77) - 5 a+hb f (u - (2a - 1)77) , (2.9) 
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while the normalization factor reads 

k-3 

f( fc )( u ) = Y[ f (u - (k - 3 - 2£)r)) . (2.10) 
1=0 

Finally, expressing T(u) from (|2.3p through the corresponding eigenvalue Q(u) one arrives to 
the formula 

T tW = Q(„ - kq) Q(u + kn) E Q(a - p /_ ^ Q(u - . - _ , (2.11, 

valid for k = 1,2,..., 00. Note, that the Bethe Ansatz equations (|2.4p guarantee that all the 
higher Tk(u) with k > 3 are entire functions of u as well as T(u). It is worth noting that 
these functions are actually eigenvalues of the "higher" transfer matrices, obtained through the 
algebraic fusion procedure [22]. In our analytic approach this information is, of course, lost. 
Nevertheless it will be useful to have in mind that the index k in the notation Tfc(it) refers to 
the dimension of the "auxiliary" space in the definition of the corresponding transfer matrix. 
Another convenient scheme of notation for higher transfer matrices (used, e.g., in [9]) is based 
on (half-)integer spin labels j, such that k = 2j + 1. 

In a generic case Eq. (]2.3p has two linear independent "Bloch wave" solutions Q±(u), defined 
by their quasi-periodicity properties, 

Q±(n + vr) = e ± *'Q±(u), (2.12) 

where the exponent ip depends on the eigenvalue T(u). These solutions satisfy the quantum 
Wronskian relation 

2*W(p)f(u) = Q + (u + T7)Q_(u-T7)-Q+(u-i7)Q_(u + 77) , (2.13) 
where W(y?) does not depend on u. Indeed, equating the two alternative expressions for T(u), 
J(u) Q+(u) = f(u - 77) Q+(n + 2t/) + f(u + 77) Q + {u - 2 V ) , (2.14) 

and 

T(tt) Q-{u) = f(u - r]) Q_(n + 2rj) + f(u + 77) Q_(u - 2r/) , (2.15) 

and writing \N(ip) as \N(ip\u) (to assume its possible ^-dependence, which cannot be ruled out 
just from the definition (|2.13j) ) one gets \N(ip\u+rj) = \N(ip\u—r]). On the other hand, Eqs. (|2.13p . 
(|2.2p and (|2.12p imply a different periodicity relation \N((p\v + w) = \N(tp\v). For generic real 77, 
these two periodicity relations can only be compatible if \N(tp\u) is independent of u, 

\N(<p\u) = \N(<p) . (2.16) 

When 77 and ip are in general position, the eigenvalues Q±(u) are locally analytic functions of 
77, therefore, by continuity, Eq. (|2.16j) at generic ip holds also when 77/^ is a rational number. 
However, when ip takes special values (for example, in the symmetric 8- vertex model) Eq. (l2,16l) 
for rational 77/^ cannot be established by the analytic arguments only. 

Obviously, the condition (|2.12p defines Q±(u) up to arbitrary u-independent normalization 
factors. Using this freedom, it is convenient to assume the normalization^] such that neither 

3 In the context of the 8 V/S OS-model this is the most natural normalization. The eigenvalues Q±(u) are 
factorized in products of theta functions (see (|3.15[) below) and the variation of ip only affects positions of zeroes. 
Obviously, the transfer matrix eigenvalues, T(u), do not have any singularities in ip. 
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of Q±(u) vanishes identically (as a function of u) or diverges at any value of (p. Then the 
quantum Wronskian W(y) will take finite values, but still can vanish at certain isolated values 
of the exponent (p. These values are called singular in the sense that there is only one quasi- 
periodic solution (|2,12p . while the second linear independent solution of (|2.3p does not possess 
the simple periodicity properties (|2.12p . As argued in [8], the singular exponents take values in 
the "dangerous" set 

Vdang = kTr + — £, k,£eZ. (2.17) 

However, each eigenvalue has its own set of singular exponents, being a subset of (|2. 1T|) . 

Evidently, Q(u) in (|2. 1 1 j) can be substituted by any of the two Bloch solutions Q±(u), so 
there are two alternative expression for each Tfc(tt). Further, multiplying (|2.14p and (|2.15|) by 
Q-(u — 2rj) and Q+(tt — 2rj) respectively, subtracting resulting equations and using ()2. 13|) one 
obtains 

2i\N(<p)T(u) = Q + (ti + 27 ? )Q„(n-2r ? )-Q + (u-2r ? )Q„(n + 2? ? ) . (2.18) 
The last result, combined with the determinant formula (|2.8p . gives 

2i\N(ip) J k (u) = Q + (u + krj) Q_(u - kr)) - Q+(u - krf) Q-{u + krj), (2.19) 

where k = 0, 1, 2, . . . , oo. 

All the functional relations presented above are general corollaries of the TQ-equation (|2.2p . 
(|231) . 



2.3 Rational values of r\ 

Let us now assume that 

2Lf] = miT, l<m<L-l, £ > 2, (2.20) 

where m and L are mutually prime integers. Evidently, 

2kr]^0 (modvr), 1 < k < L - 1 . (2.21) 

Combining the expression (|2.1ip with (|2.2p and (|2.12p . one immediately obtains the following 
functional relation, 

T L+fc (u) = 2cos(mip)T k (u + \nm) +T L _ fc (u), A; = 1,2,... . (2.22) 

which shows that, for the rational 77 of the form (|2.2U|) . all higher T^u) with k > L are the 
linear combinations of a finite number of the lower T k (u) with k < L. This relation is a simple 
corollary of the TQ-equation. It always holds for the rational values of r\ and does not require 
the existence of the second Bloch solution in (|2.12p (indeed, Eq. (|2.22p is independent of the sign 
of (p). Setting k = 1 in (|2.22j) one obtains 

Tl+i{u) = 2 cos(m ip) f(n + ^mir) + T^_i(n) . (2.23) 

This allows one to bring Eq. (|2.5p with k = L to the form 

Tl(u + 77) T L (u -v) = (% + imvr) + e imip J L -i{u)) (f(u + \nm) + e~ imip J L -i(u)) (2.24) 

where the periodicity (|2.2p of the function f (u) was taken into account. Thus, for the rational 
r/, the equations (12. 5p with k = 2,3, ... ,L — 1 together with Eq. (|2.24p form a closed system 
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of functional equations for a set of L — 1 eigenvalues |T(u), T%{u), . . . , Tl(u)}. Given that all 
T k (u) with k > 3 are recursively defined through T(u), this system of equation leads to a single 
equation involving T(it) only. Indeed, substituting the determinant formulae (12. 8|) into (12.231) 
one obtains 

det M a b(u) = , (2.25) 

l<a,b<L 

where the L by L matrix reads 

M ab {u) = M ab (u) - u S aA S bjL f(u - 3rj) - a;" 1 5 a>L 5 b)1 f(u + ??) (2.26) 
with M ab (u) given by fl^U) and to = e ±imLp . 

2.3.1 Non-zero quantum Wronskian 

Continuing the consideration of the rational case (|2.20p . let us additionally assume that both 
quasi-periodic solutions (|2.12|) exist and that their quantum Wronskian (|2.13|) is non-zero. It 
is worth noting that the functions Q±(it) in this case cannot contain complete exact strings. A 
complete exact string (or, simply, a complete string) is a ring of L zeroes u\, . . . , ul, where each 
consecutive zero differs from the previous one by 2rj, closing over the period ir, 

u k+1 = u k + 2rj, k = 1, . . . , L, u L+1 = m (mod tt) . (2.27) 

It is easy to see that any such string manifests itself as a factor in the RHS of (|2.13p . but not in 
its LHS (unless, of course, W(</?) = 0). 
It follows from (gUZD that 

Q± (u + mvr) = e ±imip Q± (u) . (2.28) 
Using (|2.19|) and (|2.20|) one easily obtains the two equivalent relations, 

e +imv T k (u)+T L _ k (u + ±mTr) = C(ip) Q + (u + kr])Q_(u-kr}) , (2.29a) 

and 

e- imi ?T k (u) + T L _ k (u+±mir) = C(<p) Q+iu-k^Q^u + kr)) , (2.29b) 
where k = 0, 1, . . . , L and 

c <"> = ■ <2 - 30) 

In particular, for k = one gets, 

T L {u + \mK) = C{ip) Q+(u)Q_(u) . (2.31) 
Quote also one simple but usefufl consequence of f|2.29j) , 

Q+ju + ky) Q+(u -krj) { e+^T k (u)+T L _ k (u+±mir) \ 

° g Q_(u + ^) 0g Q_(u-krj) ° g Tfc ( u j + j L _ k ( u + i m7r ) J [ ' > 

This first-order finite difference equation relates the ratio Q+/Q_ with the eigenvalues of the 
(higher) transfer matrices. 



4 Namely this relation with k = 1 was used in [23] to show that for rational values of r\ the expression for the 
non-linear mobility for the quantum Brownian particle in a periodic potential obtained in [8] exactly coincide 
with that of [24] found from the thermodynamic Bethe Ansatz. 
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Introduce the meromorphic functions 



uj ( „\ _ p ±im 9 V- 1 f (n + (21 + l)r?) 

* ±(t ° " 6 ^ Q ± (« + ^)Q ± ( U + (2^ + 2)^)' (2 ' 33) 

such that 

T L {u + \nm) = (Q + (n)) 2 U/+(n) = (Q-{u))' ^ '-(«) • (2.34) 

With this definition all the relations (|2,29p reduce to a single relation which again can be written 
in two equivalent forms 

nW = CM§4r' M«) = CfcO§±^. (2-35) 

Q+(«) Q-(n) 

Obviously, 

(cm) 2 , ^-(|i) 2 - C*» 

2.3.2 The RSOS regime and its vicinity 

Further reduction of the functional relation in the rational case (12.20P occurs for certain special 
values of the field from the set 

mip = (r + l)vr, r = 0,l,2,.... (2.37) 

Consider the effect of varying ip in the relation f|2.31|) . The eigenvalue Tl(u) in the LHS will 
remain finite, so as the eigenvalues Q±(u) in the RHS. The latter also do not vanish identically 
(as functions of u) at any value of <p (see the discussion of our normalization assumptions before 
(|2,17p above). Therefore the coefficient C (</?), defined in (|2.30p . is always finite. This means 
that in the rational case (I2.20p . the quantum Wronskian, W(y>), can only vanish at zeroes of the 
numerator in (|2.30p . However, the converse is not true: W(y>) does not necessarily vanish when 
C(<^) = 0. Here we are interested in this latter case where 

C(v?) = 0, W(p) / (2.38) 

with (p from the set (I2.37|) . By definition we call it the RSOS regime. The relations (|2.29p and 
(p^Tj) reduce to 

T fc (u) = (-l) r T L „ fc (n + imvr), k = 1, . . . , L - 1, (2.39a) 

and 

J L {u) = 0. (2.39b) 

All these relations can be written as a single relation (in two equivalent forms involving only 
Q_i_(n) or Q-(u) respectively), 

W + (u) = v|/_(u) = 0, (2.40) 

with v|/ ± (u) defined by (12331) . 

The special "truncation" relations (|2.39p . exactly coincide with those appearing in the RSOS- 
model [20]. These were obtained [25,26] by the algebraic fusion procedure [27] and hold for all 
eigenvalues of the RSOS model. The above analysis shows that all eigenvalues of the RSOS 
model are non-singular. The quantum Wronskian of the Bloch solutions (|2.12p is always non- 
zero (otherwise the coefficient C(ip) in (|2.29p would not have vanished). For this reason the 
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solutions of the Bethe Ansatz equations for the RSOS model cannot contain complete strings. 
Since, as argued in [13] the complete strings are necessary attributes of degenerate states, one 
arrives to a rather non-trivial statement: the spectrum of the transfer matrix in the RSOS model 
is non-degenerate. 

Consider now the vicinity of the RSOS regime, when r\ and <p are approaching their limiting 
values given by (|2.20p and (|2.37|) respectively. Interestingly, one can express some 77- and <p- 
derivatives 



d 

dr,T k {u) = — T k (u\r),<p), 



d 



calculated at the "RSOS point" , 

(77, ip) = (mir/2L, ir(r + l)/m) , 
in terms of the corresponding values of Q±(u) and their first order u-derivatives 



Q±( u ) = Q±(u\r],ip) . 



(2.41) 
(2.42) 

(2.43) 



Using (|2.19p one obtains, 



Or. 



T fe (tt) - (-iy T L - k (u + mvr/2) + Ld v T k (u) 
L 



Q' + (u + kvi) Q_(u - krf) - Q+(u - krj) Q'_(u + kr]) 



(2.44) 



J k (u)-(-lYJ L _ k (u + m7T/2) 



2\N{tp) 



Q+(u + krj) Q_(n — kn) + Q+(u — krj) Q-(u + krj) 



where the expressions in the RHS are calculated directly at the point (j2.42[) . According to the 
definitions (|2.6p . To (it) and Ti(u) do not depend on 77 and tp at all, therefore, one can express 
77- and 99-derivatives of Tx_i(u) and Tl(u) at the RSOS point (|2.42p in terms of the of values 
Q±(u) and Q±(u). 

2.4 Zero field case 

Consider now the zero field limit <p> = 0. Let us return to the case of an irrational 77 /V where 
the spectrum of the transfer matrix is non-degenerate. The eigenvalues Q±(u), corresponding 
to the same eigenstate smoothly approach the same value at cp = 0. Moreover, adjusting a 
(/j-dependent normalization of Q±(u) one can bring their small <p expansion to the form 



where 



Qo(n) = Q±(tt)|^p=o, 
From (l2~T2l it follows that 

Qo(ti + 7r) = Q (n 



Q±(«) = Qo(«) T <P QoH/2 + 0(p 2 ) 

„dQ+{u) 



Qo(«) 



dp 



dQ_(«) 



ckp 



<p=0 



Q (u + vr) = Qo(u) + 2iQ («) . 
It it easy to see that the quasi-periodic part of Qo(^) is totally determined by Qo(u), 

Q (u) = 2 -^Q (u) + Q { r\u). 

7T 



(2.45) 
(2.46) 

(2.47) 
(2.48) 
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However, the periodic part 



Qt r \u + ir) = Q { r\u), (2.49) 



can only be determined up to an additive term proportional to Qq(u). Indeed, consider the 
effect of an inessential normalization transformation 

Q ± ( u ) _, e ±ai ?Q ± (u) , (2.50) 

where a is a constant. The value of Qq(u) remains unchanged while the periodic part of Qo(u) 
transforms as 

qJTV) - Q^M - 2a Qo(u) . (2.51) 
The quantum Wronskian relation (|2.13p reduces to 

Qo(u + rj) Q (u -r,)- Q (u - rj) Q (u + r,) = 2i W(0) f (u) , (2.52) 

where 

wrm = - 

dip 



W'(0) = ^ . (2 .53) 

ip=0 



The expression ()2. 19j) now becomes 

2iW'(0)T fc («) = Q (n + kr])Q (u - krj) - Q {u - kr])Q (u + kr)) . (2.54) 

It is easy to see that at <p = the TQ-equation (|2.3p is satisfied if Q(u) there is replaced by 
either of Qo(^) or Qq(u). The same remark applies to the more general equation (|2,lip . 

The Bethe Ansatz equations (j2.4|) for the zeroes of Qo(u) are the standard equations [1] 
arising in the analysis of the symmetric 8V-model. Exactly the same equations also hold for the 
zeroes of Qo(tt), but their usefulness is very limited. Even though Qq(u) is an entire function 
of u, it lacks the simple periodicity (cf. (|2.47p ) and, therefore, does not possess any convenient 
product representation. Moreover, the transformation (|2.5ip affects the position of zeros of 
Qo(ti), making them ambiguous. All this renders the Bethe Ansatz equations for Qq(u) useless. 
Fortunately, these equations are not really required for determination of Qq(u). Once the zeros 
of Qo(fx) are known the function Qo(u) is explicitly calculated from (|2.52p (see (I3.29P below). 

Additional functional relations arise in the rational case ()2.20p . These relations are straight- 
forward corollaries of (|239jh (|23Tj) and (EOKI) . For instance, Eq. fl2T29l) gives 

T fc (u) +T L „ fc (u + 77^/2) = C(0) Qo(u + kr))Qo(u-kr)), <p = 0, (2.55) 

where 

C(0) = m/W'(0) . (2.56) 
All these relations (with different k) can be equivalently re-written as a single relation 

^Qo(n + 2^)Q (. + ( 2£ + 2 ) r ?) =C( ° ) ' ^ = °' ^ 

which is the ip = version of (|2.35p . Setting k = in (|2.55p one gets 

T L (u + mir/2) = C(0) (Q (n)) 2 , <^ = 0. (2.58) 

Thus, at ip = the eigenvalue T ^{u + mTT /2) becomes a perfect square. It only has double zeroes, 
whose positions coincide with the zeroes of Qq{u). 
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As is well known, in the rational case (|2.2U|) the transfer matrix of the 8V-model has a 
degenerate spectrum (for sufficiently large values oi N > 2L). We would like to stress here that 
the above relations (I2.55p - (|2.58j) hold only for non- degenerate states. Actually, the assumption 
made in the beginning of this subsection, that Q+(ii) coincides with Q_(u) when p = 0, is true 
only for non-degenerate states. Removing this assumption and taking ip — > limit in (|2.31j) . 
while keeping 77 fixed by (|2.20p . one obtains 

T L (u + mvr/2) = C(0) Q+(u)Q_(tt), Q± = lim Q±(u) . (2.59) 

if— >o 

For a degenerate state the eigenvalues Q+(it) and Q_(it) can only differ by positions of complete 
exact strings. This ambiguity does not affect any transfer matrix eigenvalues Tfe(tt), since the 
complete strings trivially cancel out from (|2.1ip . In principle, the complete strings can take 
arbitrary positions, however, for Q±(u) they take rather distinguished positions. Indeed, due to 
(|2.59p . the zeroes of Q±(u) manifest themselves as zeroes of Tl(u + 77i7r/2) which are uniquely 
defined even for the degenerate states. From the above discussion it is clear that T l{u) has 
either double zeroes or complete strings of zeroes. Further analysis of the degenerate case is 
contained in Sect l3.2l 



2.5 Particular models 

So far our considerations were rather general and covered several related models at the same 
time. For each particular model, one needs to specify additional properties, namely, (i) the 
explicit form of the function f(u) and (ii) detailed analytic properties of the eigenvalues Q±(u). 
In this Section we will do this for three different models: the 8V/SOS-model, the 6V-model and 
the c < 1 conformal field theory. 



2.5.1 The symmetric eight-vertex model 

The basic properties of the 8V-model are briefly reviewed in the Appendix [A] Readers who are 
not well familiar with the subject will benefit from reading this Appendix prior to the rest of 
the paper. Our notations are slightly different from those in Baxter's original papers [1—4]. The 
variables q, 77, v and p used therein (hereafter denoted as q#, t]b, vb and ps) are related to our 
variables q, 77, v and p as 

q 2 = q B = e-^, 77=^, v=™ p = p B , (2.60) 

where K# and are the complete elliptic integrals associated to the nome q#. Here we the fix 
the normalization of the Boltzmann weights as 

p=2 1^2(0 I q)" 1 ^4(0 I q 2 )" 1 , (2.61) 

where 

0i(u|q), i = 1,...,4, q = e i7rT , Imr > 0, (2.62) 

are the standard theta functions [6] with the periods ir and ttt. 

We denote the transfer matrix T and the Q-matrix from [1, 2] as T B (v) and Q B (t>), re- 
membering that our variable v is related to vb by (|2.60p . Below we often use a shifted spectral 
parameter 

u = v-ttt/2, (2.63) 
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simply connected to the variable v in (|2.6U|) . We also consider the re-defined matrices 

T(u) = (-iq- 1 / 4 )" e ivN T B {v), Q{u) = e ivN ^ 2 Q B (v) (2.64) 

where N is the number of columns of the lattice. The eigenvalues T(u) and Q(u) of these new 
matrices enjoy the following periodicity properties 

T(u + tt) = T(u), T(u + vrr) = (-q) _JV e~ 2iuN T(u) , (2.65) 

and 

8V-model: Q{u + tt) = s e inN / 2 Q(u), Q(u + 2vrr) = q~ 2JV e- 2i7Vu Q(u) . (2.66) 

Here the "quantum number" s = ±1, is the eigenvalue of the operator S, defined in (|A.12p . 
This operator always commutes with T(u) and Q(u). 

Baxter's TQ-equation (Eq.(4.2) of [1] and Eq.(87) of [2]) now takes the form with 

f(u) = (Mu\q)f ■ (2.67) 

The main reason for the above redefinitions is to bring the TQ-equation to the universal form 
(|2.3h . where T(tt) and f(it) are periodic functions of u (see Eq. (12.2p ) for an arbitrary, odd or 
even, number of sites, N. This also helps to facilitate the considerations of the scaling limit in 
our next paper [18]. 

Comparing the first equation in (|2.66|) with the periodicity of the Bloch solutions (|2.12|) one 
concludes that the exponents tp read 

(mod tt), N = even 

(2-68) 

— (mod 7r), N = odd 

Thus, for an even N the exponents of the symmetric 8V-model, with the cyclic boundary 
conditions, always belong to the "dangerous" set (|2.17p . For an odd N the exponents (|2.68p 
fall into this set only for certain rational values of 77/71". A notable example is the case r/ = tt/3, 
considered in [16-18]. 

The imaginary period relations in (I2.65P and (|2.66p certainly deserve a detailed consideration. 
First, note that in (|2,66p we only stated the periodicity with respect to the double imaginary 
period 2ttt, which always holds in all cases when the 8V-model has been exactly solveqj. Ac- 
tually, this is a rather overcautious statement which can be easily specialized further. For the 
following discussion assume a generic (i.e., irrational) value of rj/ir. Then for even N the Bloch 
solutions (I2.12|) always coincide (just as in the zero-field case of Sect l2.4l) . For odd there 
are always two linearly independent Bloch solutions for each eigenvalue T(u), one with s = +1 
and one with s = — 1 (remind that in this case each eigenvalue of the transfer matrix is double- 
degenerate [28]). The existence of the "imaginary" period imposes rather non-trivial restrictions 
on the properties of the eigenvalues. Indeed, the second relation in (I2.65P immediately implies 
that the function 

Q(n) = r q N / 2 e iuN Q(u + vrr) (2.69) 

Ref. [1] applies to rational r\ and arbitrary values of N, while ref. [2] applies to arbitrary 77 and even values 
N. It is reasonable to assume that (|2.66[) holds in general, however, the case of an arbitrary r\ and an odd N has 
never been considered. 



i 8V ) = ) 
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where r is a constant, satisfies the TQ-equation (|2.3|) as well as Q(tt). Further, if Q(u) is a Bloch 
solution 

Q(n + vr) =e* p Q(u) (2.70) 

with some ip then Q(u) is also such a solution with the exponent 

{p = ip + N7r (mod 2vr). (2.71) 

Obviously, there are two options, either Q(u) is proportional to Q(u) or it is proportional to the 
other linearly independent Bloch solution with the negated exponent "—99". The first option is 
realized for even N, 

8V-model, N even: Q(u + vrr) = r q~ N/2 e~ iuN Q(u), 

(2.72) 

The constant r = ±1 is then the eigenvalue of the spin-reversal operator 1Z (cf. (|A.37j) ). The 
second option requires the exponent (p to be a half-an-odd integer fraction of ir, it is realized for 
odd N, 

8V-model, N odd: Q±(u + vrr) = q"^ 2 e~ iuN Q T (u). 

(2.73) 

The above relations (|2.72p and (|2.T3[) were derived for irrational values of rj/ir, however they 
also hold in the rational case (|2.20p , if no additional degeneracy of the eigenvalues of the transfer 
matrix occurs (apart from the one related with the spin-reversal symmetry for odd N). The 
functional relation (12.31H can be then written in the form 

T L (-u + |mvr) = A e i7V " Q+(n) Q + (n + vrr) (2.74) 

where A is a constant. This relation is identical to the one conjectured in [14@. 

2.5.2 The solid-on-solid model 

The main idea of this paper is to study deformations of the eigenvalues T(u) and Q(u) under 
continuous variations of the exponents ip from their discrete values (|2.68[) . As explained in 
the Introduction the resulting eigenvalues correspond to the unrestricted SOS-model. We will 
therefore assume the more general periodicity relations (jl.6p for the Bloch wave solutions Q±(u), 
which hold for both odd and even N, 

SOS-model: Q±(n+vr) = e ±itp Q±(u), Q±(u+2itt) = cf 2N e~ 2iuN Q±(u), (2.75) 

where the exponent tp is arbitrary. The second exponent tp is not an independent parameter, it 
is determined by (p (see the discussion in Section [3.1.21 below). 
The second relation in (|2.75p can be further refined for even N 

SOS-model, N even: Q±(n + vrr) = q"^ 2 e^/ 2 e~ iuN Q±(u), (2.76) 

whereas the periodicity of T(u) remains the same ()2.65p as in the 8V-model. However, there is 
no a general SOS-model analog of (|2.73|) . as it is specific to half-odd exponents only. As a result 
Eq. (|2.65p is replaced with 

SOS-model, N odd: J(u + vr) = T(n), T(u + 2tit) = q~ 4N e~ 4iuN T(u) . (2.77) 

Strictly speaking the use of the term "SOS-model" here is justified for even N only [3]. Nonethe- 
less, we will use this term to indicate arbitrary values of the field parameter (p in general. 



6 The conjecture of [14] also covers a special case of degenerate states for rational values of r), where the relation 
(|2.73|) holds for the eigenvalues of the Q-matrix of [1] for even N. 
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2.5.3 Six-vertex model in a horizontal field 

The allowed vertex configurations of the six- vertex model form a subset of those shown in Figll2i 
Namely, the Boltzmann weights 07 and 08 are equal to zero. The remaining six weights will be 
parameterized as 



(2.78) 



01 = e +H - ir >a, 02 = e~ H -^a, 3 = e +H+ir ^b, 
04 = e- H+ir >b, 05 = e iu - 2ir >c, 6 = e iu ~ 2 ^c, 

where H stands for the horizontal field 

a = h(n + r ? ), b = h(u-rj), c = h(2r/), h(it) = 1 - e 2iu . (2.79) 

The above parametrization is simply related to that given in Eq.(12) of [13] (where the vertical 
field V is set to zero). The TQ-equation (eq.(ll) of [13]) takes the form ([2T2]) . ([2T31 . where 

f( u ) = (h(u)) N . (2.80) 

The Bloch solutions (|2.12p . corresponding to the eigenvectors of the transfer matrix with n 
"up-spins" , can be written as 

Q ± ( n ) = e ±W* A±(e 2m ) , (2.81) 
where A + (x) and A_(x) are polynomials in x of the degrees n and (N — n), respectively, and 

t>-*f^ + = P.82) 

Introduce new variable^] 

x = e 2iu , q = e 2lv , z = e 2ivv/n . (2.83) 

Regarding x as a new spectral parameter instead of u and writing T(u) and f(u) as T(x) and 
f(x), respectively, one can rewrite (|2.3|) in the form 

T(x)A±(x) = z ±l f{q- l x) A±(q 2 x) + z Tl f (q x) A ± (q~ 2 x) , (2.84) 

where the polynomials A±(x) are defined in (|2.8ip . This form is particularly convenient for the 
6V-model. 

2.5.4 Conformal field theory 

The continuous quantum field theory version of Baxter's commuting transfer matrices of the 
lattice theory was developed in [8,9,29]. These papers were devoted to the c < 1 conformal field 
theory (CFT). The parameters (3 and p used there define the central charge c and the Virasoro 
highest weight A, 

C = 1 - 6(/3 - /r 1 ) 2 , A= (^) +£ ^r- ( 2 - 85 ) 

They are related to our r\ and tp as 

2r] = f3 2 n, Lp = 2np/p 2 . (2.86) 



7 The parameter q should not be confused with the nome q in the 8V-model 
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The multiplicative spectral parameter A used in those papers is related to our variable u as 

A 2 = - e ~ 2iu . (2.87) 

The eigenvalues of the CFT Q-operators Q±(u) are entire functions of the variable u, satisfying 
the periodicity relation (|2.12p . Their leading asymptotics at large positive imaginary u read 

log Q±(u) = e W(T-a»») + o(i) u _> +ioo, |Re«| < vr/2 , (2.88) 

cos (^) 

where A is a known constant [8]. Here we assumed that rj does not belong to the set 

V = l(^-^r), fc = l,2,...,oo. (2.89) 



2 V 2k 

At these special values of 77 the theory contains logarithmic divergences and the asymptotics 
(I2.88|) should be replaced with 

logQ±(u) = 2i(-l) k Aue 2iuk / n + Ce 2iuk + 0(1), u^+ioo, \Reu\ < vr/2 , (2.90) 

where C is a regularization-dependent constant. The factorization formulae reacjfl 

Q ± (u) = e ±iu ^ /w A±(u), A±(«) = H ( 1 - e - 2i ("""fc) ) , (2.91) 

k=l 

where the zeroes , , . . . accumulate at infinity along the straight line 

u = ir/2 + iy, y ^ +00 . (2.92) 

Finally, the function f (u) in the case of CFT should be set to on^£|, 

f(«) = 1 . (2.93) 

With these specializations the functional relations given above become identical to those previ- 
ously obtained in [8,9,29]. 

2.6 Related developments and bibliography 

The literature on the functional relation in solvable models is huge; therefore it would not be 
practical to mention all papers in the area. Our brief review is restricted only to a subset 
of publications directly related to the eight- vertex/six- vertex models and associated models of 
quantum field theory. 

2.6.1 Transfer matrix relations 

In the above presentation the entire functions Tk(u) with k > 3 were defined by the recurrence 
relation (|2.5p . which allows one to express them solely in terms of T(u), as in (|2.8p . No other 
additional properties of Tfc(rt) were used. However, as is well known, these functions are eigen- 
values of the higher transfer matrices, usually associated with the so-called fusion procedure. 



8 Here we assumed that < r\ < 7r/4. When 7r/4 < r\ < ty/2 the product in (|2.91l) should contain the standard 
Weierstrass regularization factors [8]. 

9 Again, we have assumed that r\ does not fall into the set (|2.89[) . otherwise f(u) = exp(4 J 4?7e 2 " lfe /7r). 
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This algebraic procedure provides a derivation of the functional relations for the higher transfer 
matrices based on decomposition properties of products of representations of the affine quan- 
tum groups. Originally, all these "transfer matrix relations" were obtained essentially in this 
way. We would like to stress that the logic of these developments was exactly opposite to that 
employed in our review. The goal was to find new techniques, independent of the TQ-relation, 
rather than to deduce everything from the latter. The first important contribution was made 
by Stroganov [30] . He gave an algebraic derivation of the first nontrivial relation in (|2.5h (with 
k = 2), 

J(u + 77) T(u - rf) - f (u + 2rf) f {u - 2rf) = O ((n - u ) N ^j (2.94) 

in the vicinity of the point u = uq where the transfer matrices T(uo + rf) and T(uq — if) become 
shift operators. Remarkably, this single relation alone contains almost all information about 
the eigenvalues T(ii). To illustrate this point consider, for instance, the 6V-model. For a 
chain of the length ./V each eigenvalue T(u) is a trigonometric polynomial of the degree N, 
determined by N + 1 unknown coefficients. The mere fact that the LHS of (|2.94[) has an iV-th 
order zero immediately gives N algebraic equations for these unknowns. Similar arguments, 
obviously, apply to the 8V-model. One additional equation is usually easy to find from some 
elementary considerations (e.g., from the large u asymptotics in the 6V-model). Further, in the 
thermodynamic limit, N — > oo with u kept fixed, Eq. (|2.94p becomes a closed functional relation 
for the eigenvalues (its RHS vanishes). This is the famous "inversion relation" [30-32]. With 
additional analyticity assumptions it can be effectively used to calculate the eigenvalues of the 
transfer matrix at N = oo. Recently, Eq. (12.94p was used to derive a new non- linear integral 
equation [33], especially suited for the analysis of high-temperature properties of lattice models. 

Soon after [30] Stroganov derived [34] a particular case of (|2.23p for the 6-vertex model with 
r\ = 7r/6 (i.e., for L = 3 and m = 1). He also found that for the case of an odd number of sitea^l 
Eq. ([235|) takes the form 

J(u-2rf) T(u) T(u+2? ? ) = f (u)f {u+2rj) T(u-2r))+f(u-2r))f(u+2rj) T(u)+f (u-2ij)f (u) T(n+2r?) . 

(2.95) 

He then used this equation to obtain Bethe Ansatz type equations for the zeroes of T(u) and 
to reproduce Lieb's celebrated result [35] for the residual entropy of the two-dimensional ice. 
Unfortunately, these results were left unpublished. 

The ideas of [30, 34] were further developed in the analytic Bethe Ansatz [36] where the 
TQ-equation (or an analogous equation) is used essentially as a formal substitution to solve the 
transfer matrix functional equations. The notion of "higher" or "fused" R-matrices was devel- 
oped in [22] from the point of view of representation theory. These R-matrices were calculated 
in [37] for the 6V-model, in [38-41] for the 8V-model and in [27] for the SOS-model. The func- 
tional relations (|2,7p were given in [37] for the 6V-model and in [26] for 8V/SOS-model. The 
determinant identity (|2.25p was discussed in [26,42]. An algebraic derivation of the truncation 
relations (12.390 for the RSOS model [20] was given in [26]. A particular case of this trunca- 
tion for the hard hexagon model [43] was previously discovered in [25]. An algebraic derivation 
of (|2.23p in the zero-field six- vertex model is given in [44]. The idea of calculation of cp- and 
^-derivatives (|2.44|) at the RSOS point given in Sect l2.3.1Z1 is borrowed from [45] and [46]. 

Remarkably, the same functional equations (|2.5p (along with all their specializations in the 
rational case) arise in a related, but different context of the thermodynamic Bethe Ansatz [47]; 
see [48] for its application to the 8V-model. Usually this approach in lattice models is associated 
with non- linear integral equations. Here we refer to the functional form of these equations 

10 In our notations this corresponds to <p = ir/2 (mod n). 
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discovered in [49]. Further discussion of the correspondence of the functional relation method 
with the thermodynamic Bethe Ansatz and its generalizations for excited states can be found 
in [21,29,50-53]. 

2.6.2 Q-matrix and TQ-relations 

As noted before, a full algebraic theory of the Q-matrix in the 8V-model is not yet developed. 
The idea of the construction of the Q-matrix in terms of some special transfer matrices belongs to 
Baxter. It is a key element of his original solution of the 8V-model. Readers interested in details 
should familiarize themselves with the Appendix C of [1] (along with other four appendices and, 
of course, the main text of that paper, which contain a wealth of important information on the 
subject). The results of [1] only apply for certain rational values of ij. The construction of [1] 
and the set of allowed values of ij were recently revised in [14]. A different construction for 
the Q-matrix, which works for an arbitrary 77, was given [3]. Some remarks on a comparison 
between the two Q-matrices are given in Appendix [Al 

There are many related solvable models connected with the R-matrix of the 8V-model but 
having different L-operators and different quantum spaces. The general structure of the func- 
tional relations in all such models remains the same. In particular, they all possess a TQ-relation 
(though it may contain different scalar factors and require different analytic properties of the 
eigenvalues). In [3] Baxter also presented an extremely simple explicit formula for the matrix 
elements of the Q-matrix for the zero- field 6V-model in the sector with N/2 "up-spins" (the 
half filling). However, no such expression is known for the 8V- model, or the other sectors of 
the 6V-model. The quantum space of the 6V-model is build from the two-dimensional highest 
weight representation of U q {sli)) at every site of the lattice. Curiously enough, if this repre- 
sentation is replaced with the general cyclic representation (arising at roots of unity, q L = 1) 
then all matrix elements of the Q-matrix can be explicitly calculated [54] as a simple product 
involving only a two-spin interactional Remarkably, the resulting Q-matrix exactly coincides 
with the transfer matrix of the chiral Potts model [55-57]; this allows one to view the latter 
as a "descendant of the six- vertex model" [54]. The generalization of this construction to the 
eight-vertex and the Kashiwara-Miwa model [58] is considered in [59]. Further developments 
of the theory of the Q-matrix and related topics (along with many important applications to 
various solvable models) can be found in [60-69]. 

Baxter's original idea of the construction of Q-operators which utilizes traces of certain mon- 
odromy matrices was extended in [8,9] for trigonometric models related with the quantum affine 
algebra U q (sl(2)). It turned out that in the trigonometric case the situation is considerably 
simpler than for the 8-vertex model and the Q-operators coincide with some special transfer 
matrices. The corresponding L-operators are obtained as specializations of the universal R- 
matrix [70] to infinite-dimensional representations of the ^-oscillator algebra in the "auxiliary 
space". Although the calculations of [8,9] were specific to the continuous quantum field theory, 
the same procedure can readily be applied to lattice models (see, e.g., [71-74] for the correspond- 
ing results for the 6V- model). In the case of the 6V-model with non-zero horizontal field this 
construction leads to two Q-matriceJ^l. whose eigenvalues Q± are precisely the "Bloch wave" 
solutions of the TQ-equation. 

Note that functional relations which involves bi-linear combinations of Q±, namely (|2.19p . 

n The factorization of the matrix elements of the Q-matrix is typical for quantum space representations without 
highest and lowest weights. 

12 As noted in [9], for the "half-filled" sector of the zero-field 6V- model these Q-matrices reduce to Baxter's 
expression [3] mentioned above, as they, of course, should. 
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(|2.29|) and (|2,31|) are universal in the sense that they do not involve the model-specific function 
f(u). These relations were derived in [8,9] in the context of the conformal field theory. Similar 
relations previously appeared in the chiral Potts model [54,75,76], though the correspondence 
is not exact because there is no an additive spectral parameter in that model. Eq. (|2.19p in 
the eight- vertex and the XXX-models was found in [10] and [77]. A special case (I2.74p of the 
relation (|2.3ip involving Baxter's original Q-matrix [1] for the 8V-model was conjectured in [14]. 
Another special (zero-field) case (12.55|) of the same relation (|2.31j) in conformal field theory was 
conjectured in [52]. 



3 Quantum Wronskian relation and singular eigenvalues 

The quantum Wronskian relation (12. 131) . naturally arising in the above analysis of the TQ 
equation, is a very non-trivial functional relation. In this Section we show how this relation can 
used for the analysis of the eigenvalues. In particular, we consider a practical question of the 
calculation of Q-(u) from a known Q + (n). Next, we study certain singular eigenvalues in the 
zero-field limit. 



3.1 Solving the quantum Wronskian relation 

3.1.1 Six- vertex model (multiplicative spectral parameter) 

Let us first consider the example of the 6V-model. It is convenient to present the polynomials 
A±(x), defined in (|2,8ip in a factorized form, 

n N—n 

a+(x) = P+ n (i - x /4) > a -m = p- n ( x - x / x k) > 
fc=i k=i 

and rewrite the quantum Wronskian relation (|2,13p as 

z/K + (xq)A-{x/q) - z' 1 A+(x/g)A_ (xg) = 2i\N(z)f(x) . (3.2) 

This single functional equation determines both unknown polynomials A±(x), as well as the 
function W(z), except for an obvious freedom to choose arbitrary x-independent factors p± 
without changing the normalization-independent combination 

2i\N(z)/(p +P ^ = z-z- 1 , (3.3) 



fixed by Eq. (|3.2p at x = 0. There will be many solution to (|3.2p corresponding to different 

zt and x = g" 1 x^ 



eigenvectors of the transfer matrix. Indeed, setting there x = qxf and x = q 1 xf , one obtains 



two relations 

z A_(x+) A+{x+ q 2 ) = 2i\N(z)f(x+q), z~ x A_(x+) A+(x+ q~ 2 ) = -2i \N(z) f(x+ q- 1 ) . 

(3.4) 

Taking their ratio, one comes back to the Bethe Ansatz equations for the zeroes xf, , ■ ■ ■ , x^ 
of A+(x), 

f(x+q) 2 A + (x+g 2 ) 

z a / + _9N ' fe = l,2,...,n (3.5) 



f(x+g-i) A + (x+g-2 
Similar equations (with z replaced by z~ l ) hold for the zeroes j X2 j ■ ■ ■ ; *^7V 77, 

of A_(x) 
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Given a solution of (|3.5|) (which defines A + (x)) the corresponding polynomial A_(x) is 
uniquely determined by (|3.2p . For example, consider the case N = 2n. Comparing the leading 
power of x in (13. 2\ one gets 



S — {X l X 2 ■ ■ ■ X n ) — X± X 2 ■ ■ ■ X n . (."J-OJ 

Let us rewrite (13.21) functional difference equation 

z F(x/g) - z" 1 FW = A+i J q l X A ] +{x/q) (3-7) 

for the rational function 

. , 1 A_(x) , i 

Performing the partial fraction decomposition of the RHS of (|3.7|) and sharing the poles between 
F(xg) and F(x/q) one obtains 



F(x) 



f(x£g) (x-s 2 xl 



* V f ( x fc g" 1 ) / x-s 2 x+ \ 

(1_ S 2) A. x + A+ (x+q- 2 )A' + (x+) U-x+j 1 • j 

where A+(x) denotes the derivative of A + (x) with respect to x. The two alternative expressions 
given above are equivalent in virtue of the Bethe Ansatz equations (|3,5p . The formula (|3.9p can 
be easily generalized to arbitrary values of n. We leave this as an exercise for the reader. 

The above derivation has several limitations. Evidently the field parameter z should not 
take the values z = dbl, where the denominator in (|3.8|) vanishes. Further, the above sharing of 
the poles in (|3.7j) is only possible if all the roots x^~ are distinct, finite and non-zero. Moreover 
no pair of the roots could satisfy the condition x~j~ fx^ = q . In all these cases the expression 
(|3.9p becomes meaningless (typically the summand therein diverges). 

The quantum Wronskian relation can also be solved directly for the coefficients of the poly- 
nomials (f!T8Tj) 

n N—n 

A + (x) = ^a+x fc , A_(x)= J>,7x fc , (3.10) 

fc=0 fc=0 

This approach lead to a system of algebraic equations which sometimes easier to analyze since 
it involves only symmetric functions of the Bethe roots, but not the roots themselves. As an 
example, consider the trivial eigenvalue 

T(x) = zfixq" 1 ) + z~ 1 f(xq) (3.11) 

corresponding to the ferromagnetic ground state (all spins down) of the six-vertex. One solution 
of (|3,2p is obvious 

A+(x)=p+ = l. (3.12) 
The finite difference equation (|3.7p can be readily solved 



2»W(«) z 
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If z takes one the values z = ±g fc , k = 0, . . . , N, the polynomial A_(x) reduces to a single 
power, x k , and the corresponding Bloch solutions (|2.8ip become linearly dependent. These are 
the singular cases discussed in Sect J2.21 Indeed, choosing the normalization 

N 

P- = Y[(z 2 -q 2k ) (3.14) 
k=l 

so that A_(x) remains finite and does not vanish identically for all z and taking into account 
the definition of z in (|2.83|) . one can easily see from (13. 3p and (13.121) that the singular exponents 
form a subset of (|2.17p . as expected. 

3.1.2 Eight-vertex/SOS model 

The formula similar to (|3.9j) holds for the 8-vertex/ SOS-mo del as well. The periodicity relations 
(|2.75p imply that the entire functions Q±(u) can always be factorized as 

N 

Q±(«) =e i * ±u / ,r A ± (u), A±(u) = JJh(u-u^), (3.15) 

fc=i 

where uf,uf,..., are zeroes of Q± (u) and the elementary factors 

h(u) = t?i(u|q 2 ) (3.16) 

have the double imaginary period^!. We will assume that the roots lie in the fundamental 
domain 

0<Reu±<vr, -7t|t| < Imtifc < +7r|r|, fc = l,...,AT. (3.17) 
There are two independent sets of the Bethe Ansatz equations (12. 4p . 

with all upper and with all lower sings for the zeroes of Q+(u) and Q_(u), respectively. The 
relations (|2.75p impose the following constrains on the parameters entering (|3.15D . 

<p± = ±tp + Ntt (mod 2vr) , (3.19) 

N 

2 T tp± = =F# - 2 "J + Nir (mod 2vr) . (3.20) 

k=l 

Combining these relations one obtains 

(p + + <^__ = — 27rm,2, u+ + «_ = nmi + 2-7rrm2, mi, m.2 £ Z , (3.21) 

where 

"+ = S U k ' = U k ■ ( 3 ' 22 ) 

fc=l fc=l 

13 This choice allows a uniform treatment of even and odd N, though for even TV the zeroes Q± split into 
N/2 equidistant pairs with the separation nr. Actually, we found this redundancy quite useful in controlling the 
consistency of numerical calculations. Such pairing occurs automatically, without placing any constraints on the 
positions of zeroes. 
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Note that for any pair of eigenvalues Q±(u) the parameter ip in (|2.75|) can be regarded as a 
function of (p. Indeed, the Bethe Ansatz equations (|3.18p only contain the exponents (p±, which 
up to multiples of tt coincide with ±93. Once these equations are solved for the exponent ip 
is determined from (|3.20[) . 

With the above notations the 8V/SOS-model analog of the formula (|3.9p can be written as 

e- 2 ^+"/" h'(0) A e ^+^l- f(ut+v) h(n + 2u+-u+) 



e 



h(2u+) ^ Q + (u+ + 2r ] )Q' + (u+) h{u - u + k 



where 



h(2u+) ^ Q + «-2r?)QV«) h(«-«+) ^ J 

G( U ) = — L- (3-24) 
v ' 2i\N(ip) Q+(u) v ; 

and the prime now denotes the derivative with respect to u. This formula is a simple conse- 
quence of Bethe Ansatz equations (|3.18|) and the following elegant identity for the elliptic theta 
functional 

n $i(x-yk) _ 1 y-> I tt #i(zfc — Uj) \ fii{zk — Uk)^i(x — Zk + Z -y) 
=i - z fc ) tf^-y) fc ^ \j£ ^ (zfc _ 2i) J _ Zfc) ^ ^ 

Here a; is an independent variable, y\, . . . , y n and 21, . . . , z n are arbitrary constants, 

n n 

y = ^2vk, z = ^2z k , (3.26) 

k=l k=l 

and for brevity the theta- function #1(2; | q) is written as To obtain (|3.23p one needs to 

apply the identity (|3.25|) to the RHS of the quantum Wronskian relation, written in the form 

G( n - v )- G(« + rj) = f( "j r • (3.27) 

Q+(u + 77) Q+(ii - r/J 

The zero-field variant of (|3.23|) relates Qo(^) and Qo(u) defined in (|2.46p . Let Uk denote the 
zeroes of Qo(u), 



fc=i 

^From (|2.52p one can show that 

N 



Q (u) = Hh(u-4), ((u)= ^ (3.28) 



Qo(u) = 2 Q (n) { ™ + ^ C(« " «fc) 4(0)} + C Q (u) (3.29) 
where C is an arbitrary constant and 



TT 

k=l 



u'M = ^ = 'T k ^'\Z K r\ ■ ( 3 - 3 °) 



ifK + 77)W / (0) 
^=0 Q (n fc + 2r?)Q' (n fc ) 

The unknown constant W(0), defined in (|2.53p . can be expressed through the zeros Uk from the 
differential equations (|4.12p . 



14 In the trigonometric limit q — > this identity reduces to the modified variant of Ex.2 on page 140 of 
obtained if the constant term therein is "smartly" distributed among the terms in the sum. 
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3.2 Singular eigenvalues 

When tj/tt is a rational number the spectrum of the transfer matrix of the symmetric 8V- model is 
degenerate (for sufficiently large N) . Solutions of the Bethe Ansatz equations (|2.4p for degenerate 
states are not unique. They contains arbitrary continuous parameters, which determine positions 
of the complete strings (I2.27p . There is nothing wrong with this. In particular, this is not an 
indication that the Bethe Ansatz is "incomplete" [15]. Quite to the contrary, as explained 
in [13], the appearing continuous parameter are precisely those that are needed to describe the 
embedding of the corresponding eigenvectors into the eigenspace of the degenerate eigenvalue. 

It is not immediately clear, however, how the eigenvalues Q±(u), which have no ambiguities 
at generic values of the parameters rj and tp, could acquire continuous degrees of freedom for the 
degenerate states at special values of r] and cp. The explanation is simple: the limiting values of 
Q±(it) are not uniquely defined; they depend on the details of the limiting procedure. Here we 
will not study this phenomenon in its full generality, but give a particular example. For simplicity 
we consider the 6V-model in a field, where most calculations can be done analytically. Assume 
the same notations as in Sect J2 .5.31 In particular, recall the definitions of the multiplicative 
spectral parameter x, the field variable z and the parameter q given in (|2.83p . For generic z 
and q the transfer matrix of the 6V-model is completely non-generate (the usual degeneracy 
associated with the reversal of all spins is broken in the presence of the field). Consider the 
case N = 4 and n = 2. This sector contains six different eigenvectors. The corresponding 
polynomials A±(x), defined in (I2.81|) . are all of the second degree. Let us parameterize them as 

A + (x) = 1 + a 1 x + a 2 x 2 , A_(x) = 1 + hx + b 2 x 2 . (3.31) 

Substituting these expressions into the quantum Wronskian relation (|3.2p with p+ = p_ = 1, 
one gets four equations for the unknown coefficients a±, a 2 , &i, b 2 , 

a 2 b 2 , 

q 2 a\b 2 + Aq + a 2 b\ 
q 2 a 2 b\ + Aq + a\b 2 ' 
q% + q 2 (a 1 b 1 - 6) + a 2 

g 4 a 2 + ( ? 2 (a 1 6 1 -6)+6 2 ' [ ' 

q 2 h + 4g + ai 
q 2 ai + Aq + b\ 

Excluding a 2 ,bi,b 2 one obtains a six-degree polynomial equation for a\, which factors into two 
linear and two quadratic equations. One of the latter reads 

q(l - q 2 z 2 ) a\ + [Aq 2 (l - z 2 ) + z(l - q A )} oi + Aq{\ - z)(q 2 + z) = 0, (3.33) 

while the other is obtained by the substitution z — ► —z. Altogether one gets six different 
solutions for the coefficient a±. Once it is found, the remaining coefficients are unambiguously 
determined by (13.321) . 
Consider the limit 

rj -» tt/4, ip -> . (3.34) 

In terms of q and z it corresponds to q — ► i, z — > 1. In this limit the two eigenvalues of the 
transfer matrix, T^\x) and T( 2 )(x), corresponding to two different solutions of (|3.33p . smoothly 
approach the same value 

T(x) = 2x 4 -12x 2 + 2 = 2(x 2 -bl)(x 2 -b^ 2 ), b = 1 + V2 . (3.35) 
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To calculate Q±(u) one also needs to describe the first order deviations of r\ and ip from their 
limiting values (|3.34p . This requires two arbitrary small parameters ei and €2, or just one such 
parameter and the ratio €\/e2- We found it convenient to use the following parametrization 

v = ^ + (6 4 -6& 2 + l)e + 0(e 2 ), <p = 8(6 4 - l)e + 0(e 2 ) (3.36) 

where e — > and b is an arbitrary complex parameter, which is kept finite. The coefficients of 
the linear terms in e have been specifically chosen to simplify subsequent expressions. The cor- 
responding expansions for q and z simply follow from (]2.83[) . Substituting them into Eq. (|3.33j) . 
solving it for a\ and then determining the remaining coefficients in (|3.31j) from (|3.32j) . one 
obtains 



(3.37) 



2 

A%\x) = 1 - ^ + O(e), A W (x) = 1 - b 2 x 2 + O(e), 

2 

( z ) = 1- ^2+0(e), AL 2) (x) = l-&V + 0(e), 

6 

where 6 is related to 6 by a self-reciprocal transformation 

52 = Bw- < 338 » 

which exchanges the two solutions. The corresponding eigenvalues and eigenvectors of the 
transfer matrix read 



T«(x) = T(x) - 16(1 + 6 2 ) 2 x(l + x 2 )e + (9(e 2 



(3.39) 



T( 2 )(x) = T(x) + 32(1 - 6 2 ) 2 x(l + x 2 )e + 0(e 2 ) , 
where T(x) is given by f)3.35j) and 

= I Titi) - 1 IUt> + P (a) {| UTT} + | TTU)-| ITU)-| TUT)} + 0(e), a = 1,2, (3.40) 

with 

P 1 + 6 2 ' P 2(6 2 - 1) ' 1 ' 

Here we used self-explanatory notations for the "spin up", | |), and "spin down", | j), states of 
the edge spins. 

Note, that the case r\ = tt/4 corresponds m = 1 and L = 2 in (I2.20p . In terms of the 
variable x the complete string (|2.27p for this case consists of two roots x\,X2, constrained by 
one relation x<i = —x\. Obviously, at e = the zeroes of A±(x), given by (I3.37h . are the complete 
2-strings at x\ = b^ 1 or x\ = b^ 1 . Their positions are given in terms of the very same parameter 
b, which defines the embedding of the eigenvectors (j3.4Q[) into the eigenspace of the degenerate 
eigenvalue (|3.35p in the sector with two spins up. Remind that the parameter b can be arbitrary; 
it determines the direction of the 2- variable limit (I3.36P . 

It is interesting to see what happens in some particular cases of (|3.36p . First, consider the 
zero-field case when ip = from the very beginning and r/ is approaching the value tt/4. In 
Eq. (13.361) this corresponds b = ±1. As obvious from (|3.37p and (13.381) the limiting values of 
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A±(x) then coincide and the complete 2-string can take one of the two fixed positions with 
x\ = 1 or x\ = i. 

The second interesting case is when r\ is set to its rational value r\ = tt/4 while the ip is 
approaching zero from arbitrary values. In Eq. (|3.36|) this corresponds to b = &Q 1 , where bo is 
defined in (13.35H . Note, also the corresponding value of 6q = ^o 2 - We see that the complete 
2-strings in this case occupy precisely the distinguished positions determined by Eq. (|2.59p . In- 
deed, taking into account f)3.35|) it is easy to see that Eq. (|2.59p is satisfied with C(0) = 2 and 
Q±(a;) = A±(ir) (it does not matter which solution in (I3,37p is taken, because A^(x) = A^(x) 
when b = bo). We would like to stress that the distinguished string positions arising from (|2.59D 
can only be achieved by considering the limit from arbitrary values of the field with fixed rational 
values of r\. 

Finally, note that degenerate eigenstates in the 6V- and 8V-models were also considered 
in [13,15,78] 
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4 Analytic continuation of the eigenvalues 



The eigenvalues of the T- and Q- matrices of the 8V/SOS-model have very interesting analytic 
properties with respect to the field parameter ip. In this section we will study these properties 
by a combination of analytic and numerical techniques. We consider the disordered regime of 
the 8V-model with real q and 77 in the range 

< q < 1, vr/4 < r, < tt/2. (4.1) 

Note that the parameter r in this case is purely imaginary, r = i\r\. 



4.1 Bethe Ansatz equation 

We will show that different eigenvalues of the transfer matrix of the symmetric 8V-model can be 
obtained from each other by the analytic continuation in the variable (p. In doing this we will, 
obviously, need to consider ip as an independent variable. Therefore, for the eigenvalues Q±(u) 
we assume the SOS-model periodicity (12.75|) which allows arbitrary values of ip. In practice it 
is more convenient to work with a simply related variable, which interpolate the values of <p+ 
(remind that it is related to •p by (13.19P ). We will denote this variable by the "upright" symbol 
(p. The Weierstrass factorization for the eigenvalues was discussed in Sect J3.1.2l It is convenient 
to parameterize the zeroes of Q±(u) as 

uf = tt/2 + ia± . (4.2) 

Note, that the product representations (|3.15p has an obvious ambiguity. For example, the 
translation of any single root by the period 2tt\t\ complemented by the 2-7r-shift of tp+, 

a k^ a k~ 27 ™ r ' <p+ ^ p + + 2n , (4.3) 

leaves Q+(it) unchanged (more precisely, Q+(it) acquires an inessential scalar factor). Using this 
freedom one can always bring all zeroes of Q±(u) to the periodicity rectangular 

-tt|t| < Reaf < +7r|r|, -tt/2 < Imaf < tt/2, k = l,...,N. (4.4) 

With this convention the values of the exponents tp± in f|3. 15j) cannot, in general, be restricted 
to any finite domain. 

Introduce two functions 

* / \ 1 ! fi3{ri + ia\q) , , v 1 , $i(2r] + ia\q 2 ) 

$1 " = 7T-- ^g yj 1 r-ftr, * 2 («) = ^~ log -jj^ 1 • 45 

2tti 1/3(77 — ia\q) 2m vi(2r) — ia\q z ) 

Consider the Bethe Ansatz equations of the following general form 

]V$i(a fc ) - J] $ 3 (ait - «i) = + (4.6) 
3=1 7F 

where the numbers n\, ri2, • • • , un and (f> are given, while the (complex) variables ct\, . . . , oyy 
are considered as unknowns. The numbers take half-odd integer values for even N and integer 
values for odd N. Taking into account (|3.15p and (|4.2p it is easy to see that the logarithmic form 
of the Bethe Ansatz equations (|3.18p for Q±(u) coincide with (|4.6p provided and <j) therein 
are replaced by at and ip±. 
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Obviously, the numbers {n^} depend on the choice of branches of logarithms in (|4.5|) . In the 
considered case of the 8V-model, which involves elliptic functions, an extra care is required to 
define these branches. Since the functions are periodic 

$i(a + t7r) = $i(a), i = l,2 (4.7) 

one only needs to specify them in the periodicity strip < Ima < tt. We choose the cuts 

m + ^) 7r l r l + *(t7 + n7r ~~ v\ ( m + ^) 7r l r l + *(t7 + n7T + |> m,n£Z (4.8) 



for the function and the cuts 

{m7r|r| + i(mr — 2rf), m7r|r| + i((n — 1)tt + 2r/)}, m, n G Z (4-9) 
for $2(0) (the latter are shown in Fig. [1]). We will assume that the values of (|4.5p appearing in 
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Figure 1: The principal Ricmann sheet of the function ^(a) 



(|4.6|) are always taken on the principal sheet of the corresponding Riemann surface defined by 
the condition 

$i(-a) = -**(a), * = 1,2 . (4.10) 

The values on the cuts are taken from the right (left) side for the upper (lower) cuts in Fig(IJ as 
shown by arrows. Note that with these definitions the values of &i(a) and $2(0) on the real 
axis of a are unbounded (this does not happen in the trigonometric case). They monotonically 
decrease (or increase) with the increase of a. We should warn the reader that this property is 
never achieved with ad hoc computer definitions of the logarithms in ([4.511 . 

Every eigenvalue of the transfer matrix corresponds to some setm of the numbers {rif.}. Even 
when the branches of the logarithms are completely fixed, there is always an ambiguity in {n^} 
related to the re-definition of <p. Indeed, the transformation 

TT 2 

«fc — > a.k, n k ^n k + m, (j) ^ <p + m — , m G Z , (4.11) 



15 The numerical examples considered below suggest that these sets are the same for Q+(u) and Q_(it), provided 
the phases of the logarithms in (|4.6p are defined as described above. 
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leaves Eq. (|4.6|) unchanged. 

Further, the numbers {n k } do not uniquely determine solutions of (|4.6p . Different solutions 
corresponding to different eigenvalues could be related with the same set {n k } 0. In any 
case once the numbers n k are fixed the roots a k = o: k ((f>) solving (|4.6p become functions of 
the complex variable <f>. It turns out that the latter are multivalued functions with algebraic 
branching points. Below we will investigate their monodromy properties for various eigenvalues 
of the transfer matrix. 



4.2 Differential form of the Bethe Ansatz equations 

The behavior of the roots a k (4>) under variations of (ft can be effectively studied with a differential 
form of the Bethe Ansatz equations. Differentiating (|4.6j) with respect to (ft one immediately 
obtains a system of ordinary linear differential equations 

*Sa>-|j, (4.2) 
fc=l v 

The N by N matrix A{a)jk is given by 

N 

A{a)j k = *' 2 (ay - a*) + 5 jk {N^ l (a :i ) - ^ & 2 ( aj - a/)} (4.13) 

l=i 

where $' 12 (a) denote derivatives of the functions (|4.5|) with respect to their argument a. 
Using the identity 

flfcfr + S/lq) | G f k(x-y\g) = 6' 4 (2x\q) g 1 (2x|q)g 5 „ fc (x + y\q)6 5 _ k (x - yjq) 

fc (x + y|q) + e fc (x-y|q) 4 (2x|q) + k 2{ ™ i{ ™ d 4 (2x\q)9 k (x + y\q)6 k (x - y\q) 

(4.14) 

where k = 1,2,3,4 and 0^4 = —172,3 = 1> one can show that $' 12 (a) are meromorphic double- 
periodic function of the variable a, 

$> (g) - C v ^(Olq)^(ry + Hq^Cj - jojq) , , 

11 J ? X ^(0|q)^ 3 (r? + ia|q)^3(r/-Hq)' ^' ' 

^ + x 2^(27? + Hq 2 )^i (277 -Hq 2 )' 

where two constants £, x are given by 

, ^(2r?|q) gCOjqjgiCggjq) 

? 2vr ^ 4 (277|q) ' X 27r0 4 (2r/|q) ' 1 J 

For our purposes we need to study the differential equations (|4.12p only with the initial 
conditions corresponding to some eigenvalue of the symmetric 8V-model (typically the ground 
state eigenvalue). In this case the matrix A(a) is, in general, non-singular. Eq. (|4.12p can be then 
solved for the derivatives d,pa k ((ft) and defines locally analytic functions a k ((ft) of the variable 
(ft. Exceptions occur at certain root configurations, corresponding to special values of (ft, where 



16 With the exception of a few simple eigenvalues (e.g. the band of largest eigenvalues [79]) these sets are very 
poorly understood and their enumeration is a difficult (but, perhaps, not hopeless) problem. To our knowledge 
the only case where this enumeration problem was completely solved [80] is the c < 1 conformal field theory 
(CFT) with A — > 00. This CFT can be obtained in the scaling limit of the eight-vertex (six-vertex) model. 
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the matrix A(a) becomes singular. These singular configurations correspond to the branching 
points of the solutions in the complex <r5-plane. 
Let us rewrite (14.12j) in the form 

det\\A(a)\\^ = ^jr A %(a), (4.18) 

V k=l 

where A adi stands for the adjoint matrix and assume a(0o) = ( a i(^o)) ■ ■ ■ i <2jv(<^o)) is such that 

det||A(a(0o))|| =0. (4.19) 
The type of the branching depends on the rank of A(a((fro)). The most important case is when 

rankA(a(^ )) = N -1 . (4.20) 
Expanding the determinant det ||A(a(</>))|| near the point <j) = <j)Q, one obtains 

N 

det || J 4(a(^>))|| = ^ Q (aj(0) — aj(0o)) + higher order terms, (4-21) 
i=i 

where q denote some constants. Putting this back into (|4.18p and keeping first order terms 
only, one gets 

^{aM) ~ otiifh)) = Bj (4.22) 

i=l 

where Bj denotes the RHS of (|4.18|) evaluated at <f> = (pQ. It is easy to see that all Bj are 
non-zero, otherwise the rank of A((f>o) would have been less than N — 1. Summing (|4.22|) over j 
with the coefficients Cj one obtains 

f) N 2 N 

£T($>(«ifa) =2j2dBi. (4.23) 

^ i=i i=i 

Integrating over (j> and substituting the result back to (|4.2ip one gets 

det||A(a)|| = (2^^) ^0 - + O((0 - O )) (4.24) 

i=l 

Finally, from (|4.22p it follows that at — > 0o 

Oj-(0) - Qj(^o) = ^V</> " </>o + O((0 - &,)), £^ = 2^ (2 , (4.25) 

i=l 

Therefore, the condition (|4.20p implies the square root branching points for all otj. In principle, 
higher-order branching points are possible if the rank of the matrix A(<J)q) drops below N — 1. 
However, we found that the solutions of the Bethe Ansatz equations for N = 2,3,4 contain the 
second order branching points only. 

Note, that for even N the eigenvalues Q±(u) for non-degenerate states possess the more 
restrictive periodicity (|2.72p . Their imaginary (quasi-)period is equal to ttt rather than 2itt. 
This is also true for the SOS-model, see Eq. (l2.76p . Therefore, one can split the roots Uk{<ft) into 
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the two subsets {ai,a 2 , ■ ■ ■ ,a N / 2 } and {ajv/2) a N/2+\i ■ ■ ■ ,&n} obtained from each other by a 
uniform shift of all roots by the period tt\t\, 

a N /2 + k(<P)-a k (4>)=^\T\, k = 1,2,..., N/2. (4.26) 

It is easy to see that under this substitution the system (14.60 reduces to only N/2 equations, 
provided the phases {n N / 2+ i, • • • > n^} are properly fixed in terms of {ni, . . . , n N / 2 }. Further, the 
constraints (j4.26f) are also compatible with the differential equations (|4.12p . If these constraints 
hold for the initial data they will continue to hold for an arbitrary <f>. Therefore for even N 
one needs to consider only a half of the Bethe roots, for instance, the first N/2 roots. The 
matrix A(a) in (|4.12p effectively reduces to a N/2 x N/2 matrix. This considerably simplifies 
the numeric analysis for even values of N. 

4.3 Overview of the procedure 

Below we will give a complete classification of the eigenvalues for the symmetric 8V-model with 
N = 2, 3, 4. In most cases we present explicit analytic expressions for the eigenvalues obtained 
by a direct diagonalization of the transfer matrix. Whenever purely analytic results are not 
possible, we provide required numeric values for a particular choice of the parameters q and r\. 
Depending on convenience we will use either of the transfer matrices T(u) and T B (v) related by 
(|2.64p . assuming that the parameters u and v are always connected by (|2.63j) . The corresponding 
eigenvalues, denoted here by A(u) and A B (v) respectively, are related as 

A(«) = (-iq^ 4 ) N e iuN A B (v), v = u + ^ . (4.27) 

In the previous sections we also considered the eigenvalues of the "higher" transfer matrices 
Tfc(u), defined by (|2.5p . In this Section we will denote them as l\^ k \u). With these notations 
A(n) = A( 2 )(n). ' 

The transfer matrix of the 8V-model commutes [81,82] with the XYZ-Hamiltonian 

i N 

Kxyz = ~E ^ + J v + J* 4 J) 4 J+1) ) (4-28) 

provided the constants J x , J y , J z satisfy ()A.22p . Below we will assume the normalization 

J x = l . (4.29) 

The eigenvalues of Hxyz can be found from (|A.29P which relates this Hamiltonian to the 
logarithmic derivative of the transfer matrix T B (v). 

For each TV there are 2^ eigenvalues Aj(tt), i = 0, 1, . . . , 2 N — 1. Below we will assume a defi- 
nite orderi ng of Aj(u). We will first split them according to the momenta V of the corresponding 
eigenstatea 17 !. Then, in each sector with a fixed value of V we will arrange Aj(ii) according to 
eigenvalues of the XYZ-Hamiltonian (|4.28p . 

For every i = 0, . . . ,2^ — 1 let {a®} denote the solution of (14. 6|) . corresponding to the 
eigenvalue Aj. We will say that the contour Cij = Cij(Q, £ € [0, 1] connects two eigenvalues Aj 
and Aj (corresponding to the values fa and <f>j) if 

0^(0) = fa, C tJ (l) = fa, {«(*)} A/' ; :A r (4.30) 

17 The momenta V are eigenvalues of the shift operator, T B (r/) which cyclically shifts the edge spins a-' —* a^ 1 
to the left by one lattice step. 
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where we allowed for an ambiguity which can be compensated by an overall shift of all the 
numbers rtfc. Indeed, the only effect of the transformation (|4.1ip is that the corresponding 
eigenvalue of the transfer matrix A(tt) acquires an extra sign factor 

A(n) (-l) m A(n) . (4.31) 

As a simple illustration consider the ground state eigenvalue. In the regime (|4.ip all roots 
{ctk} are rea h so they can be arranged in an increasing sequence a± < 02 < • • • < on- For the 
ground state the numbers n\. take consecutive integer or half-an-odd integer values 

N + l 

n k = k — , k = l,...,N. (4.32) 

Further, for even N the ground state eigenvalue Aq(u) of the 8V- model is non-degenerate. It 
corresponds to <j) = 0. At this point the roots solving (|4.6p are symmetrically distributed on 
the interval (— it\t\, tt\t\), such that afc(0) = — ctN-k+i(®)- Let us take this configuration as an 
initial condition for (|4.12p and trace its evolution for real <j). When increases in the positive 
direction all the roots monotonously move in the negative direction. At <fi = ir the leftmost 
root reaches the boundary of the periodicity rectangular (I4.4h . ai(7r) = — vr|r|, while one of the 
middle roots hits the origin 0^/2+1(71") = 0. The remaining N — 2 roots arrange symmetrically 
around the origin. The value <j) = ir belongs to the set of the exponents (|2.68p of the symmetric 
8V-model. The resulting solution of the Bethe Ansatz equation corresponds to the next-to- 
leading eigenvalue of the transfer matrix, Ai(it). Thus, the contour Cq\ in this case is a straight 
line from (f) = to <f> = n. Further, at (f> = 2n the first root becomes 

ai(2vr) = a N (0) - 2vr|r| , (4.33) 

while the remaining roots take the original <fi = positions, but shifted one step left, 

a k {2ir) =a fe _i(0), k = 2,...,N. (4.34) 

The whole process is illustrated in Fig. [2j With an account of (|4.3p the resulting root configu- 
ration is completely equivalent to initial one at <p = 0. Thus the eigenvalue Aq(u) is a periodic 
function of (j) with the period 2ir. 




Figure 2: The ^-dependence of the Bethe roots for the ground state at N = 8. 

It is instructive to examine what happens to the exponent tp in (|2.75p for the same course 
variations of (f>. It can be calculated from (|3.20p . We plotted the result in Figj3l It is a 27r- 
periodic function which vanishes at the values (j> = kir, k E Z, corresponding the symmetric 
8V- model, in accordance with (|2.68|) . 
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Figure 3: The dependence between the exponents ip on <f> for the ground state eigenvalue at N = 8. 



In the remainder of this section we analyze the solutions of the Bethe Ansatz equations for 
the small size chains with N < 4. All numerical calculations are performed for the following 

case 

This choice of the parameter rj is motivated by the following considerations. First, it allows 
us to test the behavior of Bethe roots near the special point r\ = 7r/3, since the value (|4.35[> is 
quite close to it, rj — tt/3 « 0.05. Second, it provides us with all advantages of the rational case 
without running into the problems of degenerate states. Indeed, the value (|4.35p corresponds to 
m = 7 and L = 10 in (|2.20j) . while the degenerate states can only occur for N > L. Below we 
will apply the following procedure: 

(i) Analytically diagonalize the transfer matrix of the symmetric 8V-model and thus determine 
the eigenvalues Aj(it), i = 1, . . . , 2 N . 

(ii) Calculate the corresponding higher eigenvalues (it) using the formula (|2.8p . 

(iii) Use the functional relations (|2.3ip or (|2.58p with L = 10 to find zeroes of eigenvalues 
Q±(u). For even N this is straightforward since these eigenvalues coincide, however, for 
odd N one needs to correctly share zeroes between Q + (u) and Q-(u). 

(iv) Find the branching points of the eigenvalues with respect to the field variable (p. 

(v) On the Riemann surface of the largest eigenvalue Ao(w) try to find paths connecting it to 
other eigenvalues Aj(u) with the exponents from the discrete set (|2.68p of the symmetric 
8V-model. 

The case N = 2 is simple and does not really require the steps (ii)-(iii). 

The results are presented in tables of two types. In the first one we list eigenvalues of 
operators V, 5, 1Z and Hxyz (the latter are denoted as Ej) for each eigenvalue Aj(it). Remind 
that we use the normalization (|4.29p . We also present there numerical values of Kf = Af(ir/2) 
and Aj = Aj(7r/2) (remind that A B (v) and A(u) are related by (I4.27P ). The second type tables 
contain the Bethe roots and the values of the phases in (|4.6p . 

Note, that here we did not attempt to diagonalize the infinite dimensional transfer matrices 
of the SOS-model at arbitrary field and compare the results against the corresponding exact 
expressions for the eigenvalues. We hope to address this question in the future. 
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4.4 The case N 



This case is very simple and all eigenvalues can be calculated analytically. However, it is quite 
illustrative and reveals some general properties of the Bethe Ansatz equations (|4.6p . 
Four eigenvalues of the transfer-matrix Tb(v) for N = 2 are given by 

Aq = 2ab + c 2 + d 2 , Af = a 2 + b 2 + 2cd, Af = a 2 + b 2 - 2cd, Af = 2ab - c 2 - d 2 . (4.36) 

Corresponding eigenvalues of the Hamiltonian are given by 



-J X Jy ~\- J z 



J X Jy J z 



Jx + Jy + Jz- (4.37) 



For the symmetric 8V-model with even N the eigenvalues Qf{u) and (u) coincide, therefore 
there is no need to distinguish them. Their analytic expressions read 



Qo(«) 

Qi(«) 
Q 2 («) 
Qs(«) 



7T 7TT 2 7T 7TT 2 2 {0\q) 



tun i 7T , I 2\/j / 2n »2(0|g) .-tt 

e 0i(it — — + -KT\q )o\(u — —\q ) 



h(u + irT\q 2 )9i(u\q' 2 



2" ' 2 
' 2 (0|gr) _ 



j • 4 3 (it|g) 
e " 2 2 M<z) 



7TT 

~2~ 



2 0i(u|g) 



(4.38) 
(4.39) 
(4.40) 
(4.41) 



The descriptive properties of the eigenvalues are given Table [TJ Table [2] contains the Bethe 



i 


Ei 


Af 


A< 


V 




K 





-2.02064769 


1.33767166 


2.06595956 


1 


-1 


1 


1 


0.02064769 


0.66200143 


-0.04265730 


1 


1 


1 


2 


0.84245604 


0.38998228 


-0.89156871 


1 


1 


-1 


3 


1.15754396 


-0.28568795 


1.21704816 


-1 


-1 


-1 



Table 1: Descriptive properties of the eigenvalues for N = 2. 
roots, the values of the field parameter (j> = (p + = and the phases for all eigenvalues. 



i 




CU2 


4> 


m,n2 





ittt/2 


— ittt/2 





{-1/2,1/2} 


1 


ittt 





7T 


{-1/2,1/2} 


2 


mr — itt/2 


-ct/2 


7T 


{-1/2,1/2} 


3 


ittt/2 — in/2 


—ittt/2 — itt/2 





{-3/2,-1/2} 



Table 2: Roots and phases for N = 2. 

Despite being entire functions of the spectral parameter u, the eigenvalues Aj(ii) and Qi(u) 
are multivalued functions of the field variable <fi. Actually, for N = 2 their analytic properties 
are relatively simple since they are determined by the properties of only one Bethe root a\ (4>). 
Indeed, with an account of (|4.26p the Bethe Ansatz equations ()4.6p reduce to a single equation 



2$i(«i) 



2r] 



7T 



+ vr) - (m + ±) 



2/' 



(4.42) 
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which defines the function a\ (</>). Its branching points are determined by the condition (|4.19|) . 
which in this case reduces to 

=0. (4.43) 

Solving this equation for a\ and substituting the result into (|4.42p one gets potential locations 
of the branching points (f)^ in the 0-plane. 

2 

0( 6r ) = ±(f) + 2vrm + — n, m,n£Z (4.44) 
2rj 

where $' 1 (ai(<^o)) = 0. No explicit analytic expressions for <j)Q and ai(^o) are available. Their 
asymptotic expansions for small q read 



2iri sin 2r? 
7T + - 



q + q 3 (^ + ^) 



ITT /3 1 \ 

ai(4>o) = iirr iq cos 2rj + iq 3 I - cos 2rj cos 6r/ I + 0(q 5 ), q^O. (4.46) 

4 \2 6 / 

The numerical values 

O = 7T + 0.46954959i, ai(^o) = -0.727625178i + ivrr. (4.47) 

for the case (|4.35p are well approximated by the asymptotic formulae (|4.45p and ()4,46p , 

FigH] shows two sheets of the Riemann surface of the function ai(4>) corresponding to the 
ground state eigenvalue Aq(u). The principal sheet is defined by the condition ai(0) = mr/2. 
This sheet contains two cuts (cpo, it + ioo) and (cpQ, it — zoo) in the strip < Re <f> < 2-ir (they are 
shown by solid lines). The same pattern is repeated with the period 2tt. The second sheet can 
be reached from the first one by crossing either of the solid line cuts in FigU]from left to right. 
This sheet contains two new branch cuts (shown by dashed lines) in addition to the same cuts 
as on the first sheet. 

To demonstrate that all the solutions of Bethe Ansatz equations for N = 2 can be obtained 
from each other by an analytic continuation in <j) we performed numerical calculations using the 
differential equations (|4.12j) . For N = 2 it is a single equation 

dai 7] 



7r 2 <&' 1 (ai) 



(4.48) 



where &'i(a) is given by (|4.15p . 

The results are illustrated in Figj4] by three contours Coi, C±2 and C03, which connect 
corresponding eigenvalues (we used the definition (|4.30p ). The contour Coi goes between points 
on the first sheet while the contours C12 and C03 start on the first and end on the second sheet. 
Note that the ending point of the contour C03 is 7r 2 / (277) rather than as it should be according 
to the Table [2j As explained above, this does not affect the Bethe roots, but changes the sign 
of the eigenvalue as in (|4.3ip . 

A ^ -A3. (4.49) 

Apparently it is possible to find a different contour C03 ending at the correct point ^ = on 
some other sheet of the same Riemann surface. We leave this as an exercise to the reader. 

Note that in the trigonometric limit q = the vertical distance between upper and lower 
cuts vanishes and neighboring cuts join together. In the opposite low temperature limit q = 1 all 
cuts move to infinity and disappear. In both cases the Riemann surface of the eigenvalues splits 
into disjoint components, though in two different ways. In particular, for q = the connected 
parts of the Riemann surface no longer possess the 2-7r-periodicity in <p. 
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0.4695496 



-0.4695496 



4> " = S 




Figure 4: The principal sheet of the Riemann surface which contains the eigenvalue Ao for N = 2. 
4.5 The case N = 3. 

This case is quite important and will be analyzed in more details. When N is odd all eigenvalues 
of the transfer matrix of the 8V-model are doubly degenerated. So altogether there are four 
different eigenvalues. Two of them have the momentum V = 1, while for the other two V = u , 

u = e 2 ™/ 3 . 

A direct diagonalization of the transfer-matrix leads to the following analytic expressions for 
the two eigenvalues of Tb(v) with V = 1, 



A*(v) = 9i(2r,\q) 



h(v\q)Qi(v + ^o,i|q)0i(^ - ^o,i|q) 
?i(7?|q)0i(?7 + u 0l i|q)^i(77 - v ,i\q) 

where zeros Vq \ satisfy the following transcendental equations 



(4.50) 



^(0|q) 0§(«O,i|q) _ 1 + 7(2(7 " !) + 3 7 " 2C [l + 7 ± V7 2 - C7(7 + 1) + C 2 (l - 7 + 7 2 ) 



T?i(o|q)^Ki|q) 



( 7 C-l)[C(37 + 5)-3 7 -l] 



with the constants Q and 7 defined in (|A.20|) . For the choice (|4.35p its solutions rea 



7T 



v Q = 0.3749001333, «i = - + 0.29695326 i 



The remaining two eigenvalues 

0i(2»j|q)0i(t;-2»7|q)0i(v + »7|q) 2 



(v\<\) 



T lCl 



(27 7 |q) 2 6» 1 (r 7 |q)^ 1 (2^|q) 



h(v\q) 



(4.51) 



(4.52) 



(4.53) 



18 Note that Eq. (|4.51[) for 77 = 7r/3 has a solution «o = 0, corresponding to the simple ground state eigenvalue 
Aq (u) = 0i(v\q) N , considered in [16,17]. 
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Af 


Ai 


V 


s 





-1.85546785 


1.11943180 


1.40606193 


1 


±1 


1 


0.69792389 


0.25256796 


-1.45228936 


1 


±1 


2 


0.57877198 


-0.14650972 


0.65945338 




±1 


3 


0.57877198 


-0.14650972 


0.65945338 


u;- 1 


±1 



Table 3: Descriptive properties of the eigenvalues for N = 3. 
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n 1 ,n 2 ,n s 


Q+ (u) 


7T 7TT 

2 ~ ~2~ 


7T 7TT 

- + 0.42746i 

2 2 


7T 7TT 

- + — + 0.42746i 


~2 


{-1,0,1} 


Qo («) 


7T 7TT 

0.42746i 

2 2 


7T 7TT 

+ 0.42746i 

2 2 


TT 7TT 
2 + ~2~ 


TT 

+ 2 


{-1,0,1} 


Ao(u) 


_l_7rr 
~2~ 


±— - 0.374900 
2 


±^ + 0.374900 






A 3 V) 


± 

2 


±— - 1.068576 
2 


± Y + 1-068576 







Table 4: Zeros of the ground state eigenvalues for N — 3. 



correspond to T 3 = w . Their numerical values in Table [3] are real and coincide with each other 
only due to the choice of the symmetric point v = tt/2 where the second term in f|4.53j) vanishes. 
For a generic real v these eigenvalues are complex and non-degenerate. 

4.5.1 The eigenvalue Ao 

Consider the ground state eigenvalue Ao- The numerical calculation with Eq. (|2.3ip referred 
to in the step (iii) in Sect 14. 31 above require a separation of zeroes between Q + (it) and Q~(u). 
Since the total number of zeroes, 2A^ = 6, is small we used a simple trial and error method. 
Taking an arbitrary subset of 3 zeroes for Q + (n) and substituting the product (j3. 15[) into the 
TQ-equation (|2.14p one gets a relation which must be fulfilled identically in the variable u at 
some (yet undermined) value of <^+. If the subset of zeroes is chosen incorrectly then no such 
value exist and the process should be repeated. 

The numerical results (with the choice (]4.35p ) are collected in Table H] and plotted in FigEJ 
Each of the eigenvalues Q ± (u) has exactly zeros in the periodicity rectangular 

< Re{u) < tt, — k\t\ < lm(u) < tt\t\. (4.54) 

Any of the eigenvalues A^ k \u) has 2A^ zeros in the same domain (their imaginary period is ttt). 
The corresponding Bethe roots are 

.1.1 . 1 

a i — — 2 7r l r l' a 2 = 2 7r l r l — a 3 = 2 71 " ^° = 0-427465646. (4.55) 

and the field parameter cf) + = —tt/2. The roots can be obtained from by either a reflection 
ce^ = — a~^_ k or a shift 

a k = a t+i ~ n \ T \> fc = 1, . . . , AT— 1, = af + tt\t\ . (4.56) 
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Note also that the sets {a^} and {a k } precisely interlace each other (see, e.g., Fig JTUj) . The above 
remarks relating zeroes of the ground state eigenvalues Qj('u) and Qg (u) apply for arbitrary 
odd N. 
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■,tt — r\ 



TT : 



:- 



Figure 5: Zeros of Qq (u), Q (u), Aq(u), Aq 3 ^(m) for AT = 3, marked by -k , X , • , o , accordingly. 



4.5.2 The eigenvalue Ai. 

This is the smallest eigenvalue for N = 3 (it has the largest eigenvalue of the Hamiltonian) . The 
numerical zeroes (for the case (|4.35p ) are presented in Tableland plotted in Figj6j As clearly 
seen from the picture the zeroes of Qi~(it) and Q r (u) form solitary 3-strings. 
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2 


7T 7TT 
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2 2 






±Y + 1 - 1 0116874 


7TT 

T 


±— - 1.10116874 
2 



Table 5: Zeros for the smallest eigenvalue at N = 3. 

It is not difficult to check that these solutions of the Bethe Ansatz equations are connected 
with each other by an analytic continuation in <p. Consider the path from <p = 3tt/2 to <f> = —3tt/2 
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Figure 6: Zeros of (u), (u), Ai(m), A^ (w) for TV = 3, marked by ~k , X , • , o , accordingly. 



along the real axis of cp (with a small imaginary part for a better convergence). It turns out that 
the solution {a^,a~ ,a% } continued along this path transforms into another solution 

{5+} = {af - 4tt|t|,c4 + 2vr|r|,a^ + 2vr|r|}, (4.57) 

where af, i = 1,2,3 correspond to zeros of Q^(it) from the Table [5j Remind that au and Uk 
are related by (|4.2p . It is easy to check that shifts by 27rfc|r| in (|4.57p can be compensated by 
a proper choice of rik without changing the field cp + = —3tx/2. Table [5] contains the resulting 
values of n& corresponding to the roots {af , af , af}. 

The problem of the analytic continuation between different eigenvalues of the transfer matrix 
is much more complicated. The most difficult part of this procedure is a numerical study of 
the Riemann surface of the eigenvalues. For N > 3 the difficulty lies in the fact that the 
branching condition (I4.20|) just defines an (JV — l)-dimensional hyper-surface in the space of 
Bethe roots which, by itself, does not correspond to any particular value of the field (p. To 
obtain the structure of the branching points in the i^-plane one needs to determine how this 
hyper-surface intersects with required solutions of the Bethe Ansatz equations. In practice we 
used the differential equations (|4.12p and numerically studied monodromy properties of the 
solutions when (j> was varied around random loops in the complex plane. 

Remarkably, we found that the largest and smallest eigenvalues, Ao (it) and Ai (it), do indeed 
correspond to different branches of the same (multivalued) function of <j). The contour Cqi in the 
complex 0-plane connecting these eigenvalues is shown in FigJTJ As evident from the figure, the 
structure of branching cuts is extremely complicated. This figure shows the (upper half-plane) 
cuts on just two sheets of the corresponding Riemann surface. The first sheet (solid line cuts) 
contains the eigenvalue Q^(u) at 4> = vr/2, while the second one (dashed line cuts) contains the 
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eigenvalue Q^(u) at <f> = 3n/2. Note that the contour Coi leaves the first sheet when it goes 
under the cut at Re <\> = and arrives to the second of these sheets from under the (dashed line) 
cut at Re0 = tt — a. In between these two points the contour crosses additional cuts on other 
sheets of the Riemann surface, which are not shown in the picture. 




Figure 7: Structure of the cuts on the principal sheet Riemann surface containing Ao (solid lines) 
and on the sheet containing Ai (dashed lines). Only the upper half-plane is shown as the 
arrangement of the cuts is reflection-symmetric with respect to the real axis of <f> 

Note that in the limit r/ — > ir/3 the zeroes of Q^(u) form complete 3-strings which cancel 
out from the TQ-relation. As a result one gets 

Ax(u) = -6l{u + ||q) - Ol(u - ||q), r, = |. (4.58) 

The minus signs arise from the phase factors exp(±2zry^/7r) = —1 at (ft = 3n/2. When ij — > ir/3 
the parameter a tends to zero (this parameter determines positions of several cuts in FigJT]). In 
addition to that the dashed line cut with Re</> = tt and its symmetric reflection in the lower 
half plane both approach the real axes. Therefore, the contour C$\ gets pinched twice in both 
vertical and horizontal directions. As a result the complete string solution corresponding to 
Ai(u) for 7] = 7r/3 is no longer analytically connected to the ground state eigenvalue Aq(u). 

4.5.3 The eigenvalues A2,3 

It is enough to consider only one eigenvalue A2 with V = u, since all results for A3 can be 
obtained by a complex conjugation. The analytic formula for Af(v) is given in (|4.53p . The 
Bethe roots and phases for Q^(u) read 

{a^,a^ ,a^} = {-0.24023993,-0.45552631 ± 0.75825341 i}, {ni,n 2 ,n 3 } = {0,0,1} (4.59) 

and the exponent (/?_ = 7r/2. Similarly, for Qi" (u) one has 

a fc =a fc+ 7r l T l> ^+ = --y> {ni,n 2 ,n s } = {0,0, 1}. (4.60) 
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For the eigenvalue A3 all af, (p± and change their sign with respect to those for A2. 

The problem of the analytic continuation for these eigenvalues and, in particular, the question 
about their connection to the ground state eigenvalue Aq(u) have not been considered. 

4.6 The case iV = 4. 

This is the last case where we systematically study all eigenvalues of the transfer matrix. As 
before the eigenvalues will be first grouped according to their momenta, which in this case takes 
four possible values V = 1, — 1, ±i. In total there are 16 eigenvalues, all non-degenerate. Below 
we will show that all the eigenvalues with V = ±1 (there are ten such eigenvalues) correspond 
to different branches of the same multivalued function of 4> and explicitly present all paths on 
the Riemann surface, which connect these eigenvalues to each other. 

4.6.1 The sector V = 1. 

This sector contains 6 eigenvalues listed in Tabled Three of them are given by the formula 
Af = a 4 + b 4 + 2ab(a 2 + t 



+ (7 + 1)(C - 1)1 _ 2a 2 b 2 f(7 + 1)K 2 (7 " 1) " (C + l)xi] 



1-7 



( 7 -l)2 



(4.61) 

where i = 0, 3, 5 and the constants xo, X3 and X5 satisfy the same cubic equation 

x 3 + x[2C( 7 + I) 2 - (C 2 + 1)(3 - 2 7 + 3 7 2 )] + 2(C - 1) 2 (C + 1)(7 " 1) 2 (7 + 1) = 0. (4.62) 

with C, 7 defined in (|A.20|) . With the choice (|4.35p this equation has three real roots 

x = 12.98056849, x 3 = -3.41421744, x 5 = -9.56635105. (4.63) 

The remaining three eigenvalues are polynomially expressed through the Boltzmann weights 
(|A~T3|) 

Af = (a 2 + b 2 )(c 2 + d 2 ) + (a 2 + b 2 + c 2 + d 2 )(ab + cd) + Aabcd, (4.64) 

Af = (a 2 + 6 2 )(c 2 + d 2 ) + (a 2 + b 2 + c 2 + d 2 )(ab - cd) - 4abcd, (4.65) 

Af = a 4 + 6 4 -2c 2 d 2 . (4.66) 
The corresponding eigenvalues of the XYZ-Hamiltonian (|4.28p read 

E 1 = -2Ji, E 2 = -2J 2 , E 4 = -2J 3 . (4.67) 
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E< 


Af 


Ai 


V 


s 


n 





-2.75822456 


1.30674491 


2.36255547 


1 


1 


1 


1 


-2.00000000 


0.94124717 


0.45441276 


1 


-1 


1 


2 


-1.17819166 


0.61623116 


-0.92097245 


1 


-1 


-1 


3 


0.72548274 


0.14757534 


-1.33882675 


1 


1 


1 


4 


0.86310373 


0.12908441 


-1.21915423 


1 


1 


-1 


5 


2.03274182 


0.05569128 


0.61370849 


1 


1 


1 



Table 6: Propcrities of the eigenvalues with V = 1 for N = 4. 

Table [7] contains positions of the Bethe roots 01,2, the values of the field <j) and phases 
appearing in (|4.6p . Note that, since is even the roots satisfy (|4.26p and it is enough to present 
only two roots 0^2- 
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i 


«i 


Ct2 


4> 


ni,n 2 ,n 3 ,n.4 





-1.41965608 


-0.88292901 





{-3/2,-1/2,1/2,3/2} 


1 


11TT 


inr/2 


7T 


{-3/2,-1/2,1/2,3/2} 


2 


inr — in/2 


inr/2 


7T 


{-3/2,-1/2,1/2,3/2} 


3 


inr- 0.59503973 i 


0.59503973 i 





{-3/2,-1/2,1/2,3/2} 


4 


inr — in/2 








{-3/2,-1/2,1/2,3/2} 


5 


-1.65004854 -in/2 


-0.65253655 -in/2 





{-5/2,1/2,-1/2,5/2} 



Table 7: Bethe roots, values of the field and phases for V = 1 and iV = 4. 



4.6.2 The sector V = -1. 

This sector contains 4 eigenvalues. All of them are polynomial in Boltzmann weights and the 
corresponding Bethe roots can be found explicitly. 

Af = 2a 2 b 2 -c A -d 4 (4.68) 

A? = 2(a 2 b 2 + c 2 d 2 ) - (a 2 + b 2 )(c 2 + d 2 ) (4.69) 

Af = -(a 2 + b 2 )(c 2 + d 2 ) + (a 2 + b 2 + c 2 + d 2 ){ab - cd) + Aabcd (4.70) 

Af = -{a 2 + b 2 )(c 2 + d 2 ) + (a 2 + b 2 + c 2 + d 2 )(ab + cd) - 4abcd (4.71) 
The associated eigenvalues of the XYZ-Hamiltonian are 

E 6 = 2J 3 , E 7 = 0, E 8 = 2J 2 , E 9 = 2Jx. (4.72) 



i 


Ei 


Af 


A* 
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-0.86310373 


-0.51124109 


1.25718614 


-1 
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-i 
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-0.27935395 


1.63561756 


-1 
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i 
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1.17819166 


-0.09456292 


0.86905645 


-1 


-1 


-i 


9 


2.00000000 


-0.05570662 


-0.54254103 


-1 


-1 


i 



Table 8: Properties of eigenvalues for V = — 1 and JV = 4. 
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4> 


rai,n 2 ,n 3 ,n 4 
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inr/2 — in/2 


inr/2 





{-5/2,-3/2,-1/2,1/2} 
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inr/2 — i{n/2 — rj) 


inr/2 + %{n/2 - rj) 





{-3/2,-1/2,1/2,3/2} 


8 


inr/2 — in/2 





— n 


{-7/2,-1/2,-3/2,3/2} 


9 


inr/2 — in/2 


-in/2 


—n 


{-5/2,1/2,-1/2,5/2} 



Table 9: Bethe roots, values of the field and phases with V = — 1 at TV = 4. 

Using the differential equations (|4.12p we have found that all solutions of the Bethe Ansatz 
equations (|4.6p corresponding to V = ±1 can be obtained from the ground state solution by 
the analytical continuation in <j). For the numerical calculations we always assumed the values 
(I4.35p . Consider the principal sheet of the Riemann surface containing the ground state solution 
at 4> = 0. As explained before the cut structure on this sheet is always symmetric with respect to 
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the real axis and periodic with the period 2ir. Therefore we consider only the upper half of the 
periodicity strip. There is only one cut (jr + 0.67535 i, tt + i oo) on this sheet shown by solid lines 
in FigEl Encircling the branching point (as in the contour C12) brings one to the second sheet 
containing four additional cuts (shown by dashed lines in FigJS]). In this figure we also shown 
the contours which connect Ao with five other eigenvalues Ai, A2, A3, A4, — Aq. The change of 
sign for Aq is related with the transformation (|4.1ip and (|4.3ip with m = 1. The horizontal 
coordinates of the cuts shown on the figures have been numerically fitted by considering a few 
different values of r\ in the vicinity rj = 7n/2§. 




Figure 8: The principal sheet of the Riemann surface containing the eigenvalue Ao for TV = 4. 




Figure 9: The sheet of the Riemann surface containing A4 for N — 4. 



Paths to the other four eigenvalues A5, A7, As and A9 are more conveniently described on a 
different sheet of the Riemann surface which contains the eigenvalue A4 at <j) = 0. The structure 
of cuts on this sheet and connecting contours are shown in Figj9l Again at the points tt 2 / (2r]) — ir 
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and 2-7T — ir 2 /(2r]) we have to take into account (|4.31|) with m = ±1. For A7 we used the shift 
(|4.3p to move the Bethe roots to the region (|4.4p . 

4.6.3 The sector P = i. 

There are three complex eigenvalues with V = i and three complex conjugated eigenvalues with 
V = — i. Again it is sufficient to consider only V = i. 

Af = 2(a 2 b 2 - c 2 d 2 ) + i(a 2 - b 2 )(c 2 - d 2 ) (4.73) 

Afi = {ab + cd)(a 2 + b 2 - c 2 - d 2 ) + i(a 2 - b 2 )(c 2 - d 2 ) (4.74) 

Af 2 = (ab - cd){a 2 + b 2 - c 2 - d 2 ) + i(a 2 - b 2 )(c 2 - d 2 ). (4.75) 

Note that for all of them the corresponding energies Ej of the XYZ-Hamiltonian are zero. In 
general, these eigenvalues are complex. However, at the symmetric point v = tt/2 where a = b, 
all eigenvalues are real. All information about these eigenvalues is collected in Tables [lOj [11] 
The parameters of the eigenvalues with the momentum V = — i are obtained by negating of all 
ccj's, 0's and n^'s. 
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Table 10: Eigenvalues with V = i at N = 4. 
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^1,^2,713,77,4 


10 


-1.50555656 - nr/2 


-0.79702853 





{-5/2,-1/2,-1/2,3/2} 


11 


-1.89851386 - nr/2 


-1.55536378 


7T 


{-3/2,-3/2,1/2,1/2} 


12 


3fvrr/4- 0.40566532 i 


3f7rr/4 + 0.40566532 1 


7T 


{-3/2,-3/2,1/2,1/2} 



Table 11: Bethe roots, values of the field and phases for V = i and N = 4. 
The problem of analytic continuation for these eigenvalues has not been considered. 
4.7 The case N = 13 

The purpose of this section is to analyze a distribution of zeroes for the ground state eigenvalues 
for sufficiently large value of N. We deliberately choose an odd value of iV when Q+(n) and 
Q-(u) are linearly independent and their quantum Wronskian is nonzero. In this section it is 
more convenient to work with the original Baxter's normalization of eigenvalues Q±{v) and the 
variable v 

Q±(u) = e ivN/2 Q^(v), v = u + ttt/2. (4.76) 
Introduce the linear combinations Qi^iv), 

Q%(v) = (Q 1 (v)±Q 2 (v))/2, (4.77) 

such that 

Ql, 2 (« + tt) = (-1) (7V - 1)/2 Q 2 » (4.78) 
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Numerical zeroes of eigenvalues can be easily calculated from the Bethe Ansatz equations 
(|4.6p with the set of phases (|4.32p . Here we demonstrate an alternative numerical method which 
works well for the considered case. For any odd N the ground state eigenvalue A.q(v) an odd 
function of v, having N zeroes in its periodicity rectangular (|2.65p . Using the identity 

2M* + y\q)M* - y\v) = M^)Miy) -Mix)M2y), H*) = ^(||q 1/2 ) , (4.79) 

it can be represented as 

{N-l)/2 

Af(«) = 0i(«|q) E t.M^M^- 1 ' 2 ^ 2 , (4-80) 

fc=0 

where tk are constants. Similarly, Qi(v) can be written as 

(N-l)/2 

Q 1 (v)=2- N Mv) £ c k Mv) 2k Mv) N - l ~ 2k , (4.81) 

fc=0 

with some unknown constants c^. In [16] we used the representation (|4.8ip for r/ = 7r/3, however, 
it is valid for the ground state eigenvalues with arbitrary values of r\ (the case r\ = ir/3 is special 
because all coefficients can be calculated explicitly [16]). Substituting (|4.80p and (|4.8ip into 
the TQ-equation and evaluating it numerically for several values of the spectral parameter one 
gets a bi- linear system of equations for the unknown coefficients and c^. For the choice (|4.35l) 
its numerical solution is given in Table fl2l 
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5 


-18.21779332 




6 


5.231349086 



Table 12: Numerical values of the coefficients Ck in (|4.81[) for N = 13 in the case (|4.35l) . 
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-0.80200143 
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1.07830969 
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-1.30218100 


7 


-0.33007292 


11 


1.22427541 
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-1.15129255 
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0.64848502 


12 


1.38973713 




13 


1.65410007 



Table 13: Bethe roots for the ground state at N = 13. 



Using these results and (j4.76H4.78j) one can find the Bethe roots aj:, k = 1, . . . , 13. The 
results for are given in Table [TBI moreover, af = —a~^_ i+l and ip± = Ttt/2. 
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Further, zeros for Ao(u) and A { 6> {u) are as follows. The 2N = 26 zeros of Ao(u) read 

TXT 7TT 

u k = ±—, ±—±{ ri ,r 2 ,r 3) rlr^r^}, (4.82) 

where the signs can be chosen independently and r* stand for complex conjugates of rj. Similarly, 
for a[, 3) (u)> 

7TT 7TT 7T 7TT 7T 

Wfc = ±— , ±— + - ± {S1,S 2 }, ±— + - ±{s 3 ,S 4 ,S3' S 4}- ( 4 - 83 ) 

The numerical constants r j , Sj are given in Table [T4l The zeroes of eigenvalues are plotted in 
Fig.dOl 
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0.32747706 + 0.34260021 i 





Table 14: Parameters Sk, which determine positions of zeros Aq(m) and Ag (u). 
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Figure 10: Zeros of Q$(u), Q~(u), A (u), A ( 3 \u) for N = 13, marked by 



, x , • , o , accordingly. 
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Note that the sets {a^} and {a^} precisely interlace each other. This fact looks quite 
remarkable given these sets are obtained from each other by the uniform translation (|4.56p over 
the distance 7r|r| which is much larger than the spacing between roots. For a numerical study of 
the Bethe Ansatz equations on a large chain of even length (up to N = 512) we refer the reader 
to [83]. 

5 Conclusion 

In this paper we developed some new ideas in the classical subject of Baxter's celebrated eight- 
vertex and solid-on-solid models. Our primary observation concerns a (previously unnoticed) 
arbitrary field parameter in the solvable solid-on-solid model. This parameter is analogous to the 
horizontal field in the six-vertex model. This fact might not be so surprising to experts, since all 
the hard work has been done before and one just needs to lay side- by-side the papers [4, 11, 12] to 
realize that an arbitrary field parameter is, in fact, required to describe the continuous spectrum 
of the unrestricted solid-on-solid model. 

The introduction of an arbitrary field allowed us to develop a completely analytic theory 
of the functional relations in the 8V/SOS-model. The solutions of the Bethe Ansatz equations 
are multivalued functions of the field variable, having algebraic branching points. It is plausible 
that many (if not all) eigenvalues of the transfer matrix can be obtained from each other via 
analytic continuation in this variable. To demonstrate this we performed a comprehensive study 
of all eigenvalues for the 8V-model for small chains of the length N < 4 with a combination 
of analytic and numeric techniques. In particular, we saw in these cases that the largest and 
smallest eigenvalues of the transfer matrix are always connected by the analytic continuation. 
This study was partially motivated by our attempts to understand properties of the eigenvalues 
in the ferromagnetic regime [84], analytically connecting it with the disordered regime. Note 
that the ferromagnetic regime eigenvalues are important in connection with the AdS/CFT cor- 
respondence (see [85] and references therein); we hope to study them elsewhere. 

The field parameter is also important in our future considerations [18] of the quantum field 
theory limit of the 8V/SOS-model, where it becomes the massive sine-Gordon model. Note, that 
the connection between the largest and next-to-largest eigenvalues in this model was previously 
studied [53] via the thermodynamic Bethe Ansatz. The authors of [53] also considered the 
analytic continuation but in a different variable, namely, the scaling variable, which has no a 
direct analogue in the lattice theory. 

It appears that the analytic structure of eigenvalues in the eight-vertex/SOS model certainly 
deserves further studies. Somewhat simpler (but still very interesting) structure arises in the 
six- vertex model and, especially, in the c < 1 conformal field theory [80]. In the latter case the 
Riemann surface of the eigenvalues closes within each level subspace of the Virasoro module. 
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Appendix A. The eight-vertex model 



In this Appendix we briefly summarize basic properties of the symmetric eight-vertex (8V) 
model used in this paper. For a more detailed information the reader should consult with 
Baxter's original publications [1—5]. 

Consider a square lattice of N columns and M rows, with the toroidal boundary conditions. 
Each edge of the lattice carries a spin variable a taking two values a = + and a = — , corre- 
sponding to the "spin-up" and "spin-down" states of the edge. Each vertex is assigned with 
a Boltzmann weight R(a, a'\(3,(3') depending on the spin states a,a',/3,/3' of the surrounding 
edges arranged as in Figllll There are only eight "allowed" vertex configurations, shown in 



R(a,a'\/3,/3') = 







a 


a' 







Figure 11: The arrangement of four spins around the vertex. 

FigH21 which have non-vanishing Boltzmann weights. These weights are not arbitrary; they 
parameterized by only four arbitrary constants a, b, c, d, 

01 = 02 = a, 03 = 04 = b, 05 = 6 = C, 7 = 8 = d . (A.l) 

The remaining eight configurations are forbidden; their Boltzmann weight is zero. 
The partition function 

Z = E II R(<*,°{\Prf) (A.2) 

(spins) (vertices) 

is defined as the sum over all spin configurations of the whole lattice, where each configuration 
is counted with the weight equal to the product of vertex weights over all lattice vertices. 

The vertex weight matrix R(a,a'\f3, f3') can be thought as a matrix acting in the direct 
product of the two two-dimensional vector spaces C 2 (g) C 2 , where the indices a, a' refer to the 
first space and the indices (5, (5' to the second. It has an elegant representation [1] 

4 
k=l 

where cr^, k = 1,2,3,4, are the Pauli matrices 




(A.4) 
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Figure 12: Eight allowed vertex configuration and their Boltzmann weights. Thin lines represent 
the "spin-up" states and the bold lines represent the "spin-down" states of the edge 
spins 



and 

c + d c — d a — b a + b / A ^\ 

wi = — ^— , w 2 = — ^— , w 3 = — ^— , w A = — ^— . (A.5) 

The matrix R can conveniently be presented as a two-by-two block matrix acting in the first 
space with the two-dimensional matrix blocks acting in the second space, 

where 

The matrix R possesses simple symmetry relations, obvious from (1A.3[) . 

R = (a k ®a k ) H(a k (g)a k ), k = 1,2,3. (A.8) 

The row-to-row transfer matrix T acts in space of states of iV spins located on a horizontal 
row of vertical edges. It is defined as the trace of the product of the two-by-two matrices (|A.6|) 

T = Tr c2 fRWR^-.RW), (A.9) 



where the operator entries (IA.7j) of the matrix R") act on the states of the j-th spin in the row. 
The partition function (|A.2[) can written as 

Z = Tr(T M ) , (A.10) 

where M is the number of rows of the lattice. It follows from (|A.8[) that 

[T,S] = [T,K] =0 (A.ll) 



49 



where the operators 



S = G 



K = a 



(i) 



a 



(2) 



(A.12) 



(i) (i) 

defined as the product of the local spin operators and cr} over the all sites, j = 1, . . . , N , 
in a row. Note, that for an odd iV these operators do not commute among themselves (indeed, 
TZS = (-1) N SK) and only one of them can be diagonalized simultaneously with T(v). Below 
we will always assume the basis where the operator S is diagonal (then for even TV the operator 
1Z will be diagonal as well). 

Following [1] we parameterize the Boltzmann weights a, b, c, d as 

a = p -#4(277 j q 2 ) #4(7; — 77 1 q 2 ) #i(u + 77 1 q 2 ), 

b = p #4(277 I q 2 ) -&i(v - 77 1 q 2 ) #4(7; + 77 1 q 2 ), 

c = p #1(277 I q 2 ) #4(7; — 77 1 q 2 ) #4(1; + 77 1 q 2 ), 

d = p #1(277 1 q 2 ) #1(7; — 77 1 q 2 ) #1(7; + 77 1 q 2 ), 



(A.13) 



where we use the standard theta-functions [6] 

$i(v\q), i = 1,...,4, 



Imr > 0, 



(A.14) 



with the periods ir and ttt. Apart from the simple change in the notation for the theta-functions 
the above parametrization is the same as that given by Eq.(8) of [1]. The theta-functions H(z) 
and Q(z) of the nome qs therein are given by 



q 2 , H(z)=0i( 



7TZ 



q 2 ), Q(z)=M 



irz 



2K B 1 " " w y 2K B 1 ^ ' 
where K# is the complete elliptic integral of the first kind of the modulus 

#2(0 I q B ) 



k 



# 3 (0 I q B ) 



(A.15) 



(A.16) 



The variables v, 77 and p used in [1] (hereafter denoted as vb, r]B an d Pb) are related to those 
in (LAT3I) by 

irv B irr]B , A 17 , 

U= 2fV V= 2Kb-> P = PB - (A - 17) 

With these definitions the transfer matrix (|A.9P is the same as that given by Eqs.(6)-(8) of ref [1]. 
We will denote it as T B {v), remembering that our variable v is related to Baxter's variable vb 
by fATTl) . 

Throughout this paper we fix the normalization factor p in (|A.13[) as 



p = 2 7f 2 (0 I q)- 1 ^ 4 (0 I q 2 )" 1 . 
Note, that the variables (|A.5P can be then written as 



Wi = -#l(2?7 q) p — j — r, 7 

2 # 5 _i(77|q) 
We shall also make use of two invariants 



1, 



,4. 



c 



cd 
oh 



#i(277|q s 
#4(277|q 2 



7 



^3 ~ w 2 

w\ — w\ 



fli(T/|q)fl 4 (r;|q) 
e 2 (77|q)0 3 (r/|q) 



(A.18) 



(A.19) 



(A.20) 
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The transfer matrix T B (v) commutes [81, 82] with the Hamiltonian of the XYZ-model 
(remind that we are assuming the periodic boundary condition), 



1 N 

Kxyz = ~E V* + Jy 4 3) 4 3+1) + J* 4 j) <rt 1] ) (A. 21) 



j'=i 
provided 



^ 4 (2t/ | q) . ^ 3 (2r7 1 q) _ # 2 (2r7 1 q) 
i? 4 (0|q) ' i? 3 (0|q) ' i? 2 (0|q) 



J*:Jy:J,= ^VZV ■ ^VZV ■ ^TTT ■ (A.22) 



Following [82] let us show that this Hamiltonian is simply related to the logarithmic derivative 
of the transfer matrix at v = rj. It follows from (|A.3|) and (|A.19|) that 



, AT 

^ 1 (2r ? |q) f logT» = iVp' 4 I+£ a««7? +1 )+^ a?V? +1 >+p& ^Vf +1) ) , (A.23) 

LOU J v=rt * — ' 



where I denotes the unit operator. Here we used new variables 

pi = - w 2 - w 3 + W4), P2 = ^(~W! + w 2 - w 3 + w 4 ) , (A.24) 

P3 = ^(-wi-w 2 + w 3 + w 4: ), p A = -{wx + w 2 + w 3 + W4) . (A.25) 
Their ^-derivatives evaluated at v = r/ denoted as 



, dpi 

p *= *7 



v=t) 



(A.26) 



Note that 

Pi +P2 +P3 +Pa = Zwi . (A.27) 
Using (|A.19j) on can readily show that 

2pi = ^(0\ q / 5 - i{ ^ , , = 1,2,3. (A.28) 
^5-i(0|q) 

Combining the last relation with (|A.23|) one can express the Hamiltonian (|A.21j) as 

H-xyz = - .^ff^l S {M2v I q) £logT» - Np' 4 l\ (A.29) 
174(277 1 q) i?i(0 1 q) I dv ) v=r) 

In the main text of the paper we mostly use the spectral parameter u which is related to 
variable v in (|A.13P as 

u = v- vrr/2 . (A.30) 

In (|2.64p we also define the renormalized transfer-matrix T(u) differing from T B (u + ttt/2) by 
a simple u-dependent factor. Of course, this factor could have been included in the normaliza- 
tion of the Boltzmann weights (|A.13p . but we prefer not to do so and keep a clear distinction 
between two alternative normalizations and, in fact, to use both of them. In particular, the 
Baxter's normalization and the variable v are more suitable in writing explicit expressions for 
the eigenvalues of the transfer matrix presented in Sect. [H 
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There are two related, but different, constructions of the Q-matrix for the 8V-model, given 
in Baxter's 1972 and 1973 papers [1] and [2]. The results of [1] apply for rational rj's and an 
arbitrary number of sites N, while the those of [2] apply to arbitrary 77's and even N. We will 
denote the corresponding Q-matrices as Q^ 72 ) (v) and Q^ 73 ) (v). They obey the same TQ-equation 
(Eq.(4.2) of [1] and Eq.(87) of [2]), 

T B (v) Q B (v) = </>(v - 77) Q B (v + 2r?) + <f>(v + 77) Q B (v - 2rj), (A.31) 

where 

cf ) B (v) = (Mv\q)) N • (A.32) 
and possess the same periodicity properties 

Q B (v + ir) = SQ B (v), Q B (v + 2irT) = q- 2N e- 2lvN Q B {v). (A.33) 

where S is defined in (jA. 12|) and Q B (v) stands for either of Q^ 72 ^^) or Q( 73 )(i;). Nevertheless, 
as noticed in [14], these two Q-matrices do not coincide. Of course, they can only be compared 
for rational 77 and even N when both of them can be constructed. As shown in [14] the matrices 
Q*- 72 )(?j) and Q( 73 )(tj) have different eigenvalues for degenerate eigenstates of the transfer ma- 
trix. This phenomenon is obviously related with the non-uniqueness of the eigenvectors for the 
degenerate eigenstates [13] and demonstrates the impossibility of an universal definition of the 
Q-matrix in the zero- field 8V-model at rational 77 (cf. [54]). 

The difference between Q( 72 )(?j) and Q( 73 )(tj) can be traced back [14] to their commutativity 
properties with the operator 1Z. The matrix Q( 73 )(t;) (defined for even N only) always com- 
mutes with 1Z. The matrix Q( 72 )(t;) (defined for all values of N) does not commute with 1Z, 
irrespectively to whether N is odd or even. On the other hand both these matrices commute 
with S, so that the first equation immediately translates into the (real period) periodicity for 
the eigenvalues (the first relation in (|2.66p ). The second equation in (|A.33|) implies only the 
2-7rT-periodicity relation (the second relation in (|2.66j) ) . 

To complete our discussion of various constructions for the Q-matrix, note that in [5] Baxter 
considered a modified version of his Q^ 73 ^ matrix. There he used yet another set of parameters, 
namely the variables v and A, which we denote here as ' and A( 82 ). They are related to v 
and rj in (|A.13j) as 

v (82) = ^B (2v _ v) A (82) = 2iKB (?r _ ^ (A 34) 

7T 7T 

where Kb is the same as in (fATTTl) . Writing the Q-matrix of [5] as Q( 82 )(t/ 82 ), A (82) ) and that 
of [2] as Q( 73 )(w,r?), and using the explanations in Sect. 7 of [13] it is not difficult to verify that 

Q( 73S >(v,7]) =KPQ( 82 >(ii( 82 ',A (82) ) (A.35) 

where the operator 

V = e +i r fJ3 (8) e _i f CT s ® e + ^ as ® • • • (g> e~^ as (A.36) 
is defined similarly to those in (|A.12|) . Using the relations 

nV = {-l) N ' 2 VlZS, QW(v + TrT)=q- N / 2 e- ivN KQ( 73 \v), (A.37) 

one can show that 

Q( 82 V 82 ) +2K' b ) = (q B r m e W ^v m N/AK B )nSQ^ 2 \v^) (A.38) 
which is exactly Eq. (10.5.43a) of [5]. 
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The partition function (|A.2|) can be regarded as a function of four parameters w\ , w 2 , w 3 , W4 , 
defined in COl) . S o one can write it as Z{w\, W2, w 3 , W4). When both M and N are even, it 
possesses the symmetry 

Z(wi, W2, w 3 ,wa) = Z(±Wi,±Wj,±w k ,±wi), (A.39) 

where {i,j, k,l} is an arbitrary permutation of {1,2,3,4} and all signs can be chosen indepen- 
dently. Using this symmetry one can always rearrange Wi so that 

W\ > W2 > W 3 > 1 104 1 . (A. 40) 

The partition function-per-site, k(v), is defined as 

log k(v) = A/ lnn ^ ^ log Z = I log A » (A.41) 

where A B (v) is the maximum eigenvalue of the transfer- matrix T B (v). In the (unphysical) 
regime (TA40]) it reads [1] 

log K (v) = log c + 2 V ^h 2 ((f B - f B )n) (cosh(nZ 5 ) - cosh( TOij )) 

n sinh( 2nfR ) cosh ( nZ r ) 

n=l 

where 



as = —2iv, Ib = —2ir], tb = — hit. (A. 43) 

In this paper we consider the disordered regimj^l 

< rj < 7r/2, 7]<v<ir-r], (A. 44) 

which corresponds to a different ordering of the variables Wi, namely, 

Wi > wi > w 2 > w 3 > 0. (A.45) 

This regime can be mapped into (|A.40|) by the following transformation of the weights 

, a — b + c — d ,, a — b — c + d . a + b + c + d a + b — c — d , . , „. 
a = , b' = , c' = , d' = , (A.46) 

which is equivalent to 

w[ = io 4 , w' 2 = wi, w 3 = w 2 , w' A = w 3 , (A.47) 

where idJ,^, w' 3 , w' a are defined by (|A.5|) with a, b, c, d replaced by a', b' , d , d' . It is easy to check 
that if wi , W2 , w 3 , w 4 belong to the regime (|A.45j) then w[ , w' 2 , w' 3 , w'^ belong to (IA.40j) . Putting 
(TAT31) into (|AT46l) it is not very difficult to show that the "transformed" weights a',b',c', d! can 
be parameterized by same formulae (IA.13jl provided one makes the following substitution 

V-TT/2 77 - 7T/2 1 .. A . 

r r r 

where 

pi = -i(-ir)" 1/2 exp{ — —[v(v - w) + 377(77 - ir) + ir 2 ] f . (A.49) 



19 See [5] for the classification of the regimes of the 8V-model. Beware that the variables therein are numbered 
differently as compared to this work. Here we follow the original papers [1-4]. 
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Thus, to find the partition function-per-site in the regime (|A.45p one needs to substitute pa- 
rameters ob,Ib an d tb in (|A.42|) by 

a B = , Ib = , T B = — , (A.50) 

T T T 

and replace c by d . After straightforward calculations one obtains 

2in(v — n) xr^ sinh(2r7x n ) sinh(2(v — 77) x n ) ,, 

\og K (v )=log& 1 (v + r, q + — ~ + E ■ u \ UT TY^ ' ( A - 51 ) 

7rr n smh(7rx n J cosh((7r — In) x n ) 

where x n = injr. The above formula can be equivalently rewritten as 

ir](ir-3rj) ^ sinh(r/ x n ) sinh((7r - 3??) x n ) cosh((7r - 2«) x n ) 

log«(u) = logMv q) + + 2 > — r TTT ^ v • 

7rr n smh(7rx n J cosh((7r — 2,-q) x n ) 

(A.52) 

Note that for i] = n/3 the second and third terms in the last expression vanish, so it simplifies 

to 

log k(v) = log "&\{v I q), r] = 7r/3. (A. 53) 

Recall that the above derivation was based on the symmetry (|A.39P valid for even values of N 
only. However, since in the disordered regime (IA.45j) the Boltzmann weights (1A.13D are strictly 
positive, the result (|A.53p . obviously, does not depend on the way the limit N — > 00 is taken. 
Eqs. (|A.41|) and (IA.53P imply that the asymptotics on the largest eigenvalue of T B (v) in this 
special case has an extremely simple form Aq(v) ~ ($i(v \ q))^. Remarkably, as it was argued 
in [16,42,86], for odd values of this asymptotics is exact even for a finite lattice, 

A^(v) = {a + b) N = ^(v\q), N = 2n + 1, 77 = tt/3 . (A.54) 

We have also checked this fact analytically up to = 15. 
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